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Abstract: We study the thermal partition function of level k U(N) Chern-Simons theories 
on S 2 interacting with matter in the fundamental representation. We work in the 't Hooft 
limit, N, k — > oo, with A = N/k and T ^ 2 held fixed where T is the temperature and 
V2 the volume of the sphere. An effective action proposed in arXiv: 121 1.4843 relates the 
partition function to the expectation value of a 'potential' function of the S 1 holonomy in 
pure Chern-Simons theory; in several examples we compute the holonomy potential as a 
function of A. We use level rank duality of pure Chern-Simons theory to demonstrate the 
equality of thermal partition functions of previously conjectured dual pairs of theories as 
a function of the temperature. We reduce the partition function to a matrix integral over 
holonomies. The summation over flux sectors quantizes the eigenvalues of this matrix in 
units of ^ an d the eigenvalue density of the holonomy matrix is bounded from above by 
2^r. The corresponding matrix integrals generically undergo two phase transitions as a 
function of temperature. For several Chern-Simons matter theories we are able to exactly 
solve the relevant matrix models in the low temperature phase, and determine the phase 
transition temperature as a function of A. At low temperatures our partition function 
smoothly matches onto the iV and A independent free energy of a gas of non renormalized 
mufti trace operators. We also find an exact solution to a simple toy matrix model; the 
large N Gross- Witten-Wadia matrix integral subject to an upper bound on eigenvalue 
density. 
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1. Introduction 



The AdS / CFT correspondence maps deconfinement transitions of large N gauge theories 
on spheres to gravitational phase transitions involving black hole nucleation Efl. This 
observation has motivated the intensive study of deconfinement phase transitions of large 
N p+1 dimensional Yang Mills theories (coupled to adjoint and fundamental matter) on 
S p . The thermal partition function of a Yang-Mills theory on S p is given by the Euclidean 
path integral of the theory on x 5 1 . Upon integrating out all massive modes this path 
integral reduces to an integral over the single unitary matrix U 



N 

Z YM = / DU exp[-V YM (U)} =11 

m=l ~ 



da r 



Y[ 2 sin | : | r 



a l~ a m \ -V YM (U) 



to, 



where U is the zero mode (on S p ) of the holonomy around the thermal circle, e 
1 . . . N) are the eigenvalues of U. Vym(U) is a potential function whose precise form 
depends on the theory under study. At least in perturbation theory [|| and perhaps beyond 
g, (§, the potential V YM (U) is an analytic function of U. (|1,1D may be thought of as 



a Landau Ginzburg or Wilsonian description of the holonomy U, the lightest degree of 
freedom of the finite temperature field theory. 

The effective potential Vym{U) was computed in free gauge theories ||, ^]; it has also 
been evaluated at higher orders in perturbation theory in special examples || 0, g, ||. 
In all these examples Vym(U) is an attractive potential for the eigenvalues of the unitary 
matrix, t' Hooft counting and the requirement of extensivity force Vym(U) to scale like 
N 2 V p . In the special case of a conformal theory the potential scales like N 2 V P T P and the 
attraction between eigenvalues grows arbitrarily large at high temperatures. 1 On the other 
hand the integration measure DU (which vanishes when any two eigenvalues coincide) 
supplies a temperature independent repulsive potential for the eigenvalues. At large N 
the leading piece of the partition function is determined by a saddle point distribution of 
eigenvalues of U. Repulsion from the measure dominates over attraction from the potential 
at low temperatures and the eigenvalues of U are distributed all over the unit circle in the 
complex plane. 2 The integral ( |1.1[ ) undergoes Gross- Witten-Wadia type deconfinement 
transitions |i~0| , 11, 12] at V p T p of order unity. 3 At high enough temperatures the saddle 



point eigenvalue distribution of U has support on only a small arc in the unit circle of the 
complex plane. The size of this arc goes to zero as V p T p — > oo. In this 'decompactification' 
limit, the holonomy matrix U is localized around the unit matrix. 

In this note we study the finite temperature phase structure of renormalized level k 
U(iV) Chern-Simons theories coupled to a finite number of fundamental fields on S 2 in the 

1 Most of the general discussion of this paragraph applies also to non conformal theories in the appropriate 
high temperature and/or large volume limit. 

2 If the matter content of the theory consists only of fundamental plus adjoint fields, and if the matter 
content is held fixed as N which is taken to infinity and the low temperature saddle point for U is the clock 
matrix; a matrix whose eigenvalues are uniformly distributed on the unit circle in the complex plane. 

3 At weak coupling at least the system undergoes either a single first order transition or a second order 
transition followed by a third order phase transition depending on the details of a quartic term in the 
potential V{U) §. 
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t' Hooft limit N — > oo, k — > oo with A = fixed. 4 We address the following question: Does 
the Landau Ginzburg description ( |1.1| ) continue to hold for Chern-Simons matter theories? 
If not what replaces it? Our analysis of this structural question applies to all fundamental 
matter Chern-Simons theories in the high temperature and large N limit; in addition we 
make significant progress towards an exact determination of the partition functions the 
large N free energy, as a function of temperature and A, in particular examples of these 
theories. 

To start our discussion let us first review the behaviour of matter Chern-Simons the- 
ories in the free limit A — > 0. In particular examples such a thermal partition function 
was studied in |l6|, [l?]]. The free partition function continues to be governed by the 
Landau Ginzburg form ( |l.l|) even at finite N; in fact the only qualitative difference from 
the Yang Mills case is that V{U) that scales like NV 2 T 2 (rather than N 2 V 2 T 2 ) at high 
temperatures. 5 As the repulsion from the measure DU in ( |1.1| ) continues to scale like N 2 , 
the partition function displays interesting dynamics in the large limit only if we set 
V2T 2 = £iV and take the limit N — > 00 holding ( fixed. At small ( the eigenvalue density 
of the holonomy matrix has support everywhere on the unit circle in the complex plane; 
however at £ of order unity the system undergoes a third order clumping or deconfinement 
phase transition. At higher temperature the eigenvalue distribution has support only on a 
finite arc (centered about unity) on the complex plane. In the infinite £ limit the size of 
this arc goes to zero and the holonomy reduces to the identity matrix. 

Large N fundamental matter Chern-Simons theories have recently been studied inten- 



sively at finite values of the t' Hooft coupling A (see e.g. |16L Eq, [jjj, B0[ ^T], [L7|, [22, £3 



24, 25, 26, |27|]). In particular it has been argued in the important recent paper [27] that 
the entire effect of matter loops on gauge dynamics, at temperatures of order y/~N and at 
leading order in N, is to generate a term of the form 



T I d 2 x s/gv{U{x)) 

for the gauge theory effective action. Here v (U) is an effective potential of order iV, 6 whose 
detailed form depends on the matter content and couplings of the Chern-Simons matter 
system. It follows that in the large N limit, the thermal Chern-Simons matter path integral 



4 We use the dimensional reduction regulation scheme through this paper. In this case |fcj = \kym\ + N, 
where kyu is the level of the Chern-Simons theory regulated by including an infinitesimal Yang Mills term 
in the action, and k is the level of the theory regulated in the dimensional reduction scheme ^| . For 
this reason this class of Chern-Simons theories may be well-defined only in the range |A| < 1. Hereafter we 
assume the 't Hooft coupling is always positive for simplicity. Our results are easily generalized to the case 
A is negative by taking the absolute value of it. 

5 The potential scales like N rather than iV 2 because the matter is in the fundamental representation 
and pure Chern-Simons theory has no propagating degrees of freedom. 

6 The trace of U is counted as order N. An example of a potential of order iV is 

v(U) = aTrU + b ( ' +... + c.c. 

N- 1 
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is given by 



Zcs = / DA exp 



d 2 xj~g v(U(x)) 



(1.2) 



( |1.2| ) is simply the pure Chern-Simons theory with an added potential v(U(x)) the local 
value of the holonomy matrix U{x). As pure Chern-Simons theory is topological, the 
expectation values of observables are independent of x and (1.2) may be rewritten as 



Z CS = (e- T2v ^) N , k (1.3) 

where {0)^,k denotes the expectation value of O in the pure Chern-Simons theory at rank 
iV and level k and V2 is the proper volume of the spatial manifold (S 2 in this case). In 
other words the partition function of the matter Chern-Simons theory may be re expressed 
as linear combination of Wilson loops (in various representations of U(N)) that wind the 
time circle. The only memory of the fundamental matter lies in the form of the function 
v(U), whose structure is determined by the fundamental matter content and interactions. 

Before turning to the important question of how v (U) may be determined in any given 
matter Chern-Simons theory, we pause to note an important property of the formula (|L3j). 
As Ql.3| ) is the expectation value of a (sum of) Wilson loops in rank N and (renormalized) 
level k Chern-Simons theory, the well established level rank duality of pure Chern-Simons 
theory relates the expectation value in (|l.3| ) to the expectation value of a dual operator 
in rank k — N and level k pure Chern-Simons theory. 7 (see p8| , ^?J] for a relatively recent 
discussion of level rank duality and references to earlier work). The specific relationship 
turns out to be the following. Any gauge invariant function v (U) may be regarded as a 
function of the N variables TrU n . For any such function we define a corresponding dual 
function v(U) by the equation 

v {trU n ) = v ({-l) n+1 trU n ) . (1.4) 
Level rank duality of pure Chern-Simons theory turns out to imply that 

\e ')N,k — \e ')k-N,k- (,1-oj 



Ql.5| ) will have interesting implications for matter Chern-Simons theories as we will see 
below. 

Let us now return to computation of v(U) for any given Chern-Simons matter theory. 
It is a remarkable fact that the function v(U) appears to be exactly computable as a 
function of A in the large N limit for arbitrary matter Chern-Simons theories. As v (U) is 
a local expression, it may be computed on R 2 x S 1 rather than S 2 x S , Using standard 



large N techniques and employing an unusual lightcone gauge, the authors of [16| were 
able to compute v(U) for the special case that U = I as an arbitrary function of A for 
the theory of fundamental fermions minimally coupled to the Chern-Simons field. This 



7 If we use the unrenormalized Chern-Simons level kyM then the level rank duality states the invariance 
of pure Chern-Simons theory under the exchange of kyM and TV. 
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computation was later generalized to other theories in [23|,|}(]]. It was also generalized to 
the computation of v[U] for a matrix U whose eigenvalue distribution is given by ( 1.1C| ) 
below (more below on why this is important) in p7[| . A straightforward generalization of 
these computations permits the evaluation of v[U] for arbitrary U; our results for v[U, A] 
are presented in section |3| below for several examples of matter Chern-Simons theories. 
Pairs of the matter Chern Simons theories we have studied in section |3| below have been 
conjectured to be related via level rank type dualities. In section ^ below we demonstrate 
that, in each case, the potentials v[U] for conjecturally dual Chern-Simons matter theories 
are related by fll.4|) . It follows from (|l.5|) above, that S 2 partition functions for conjectural 
dual matter Chern-Simons theory pairs are equal at every value of A and temperature. 



Such dual pairs include Giveon-Kutasov duals for a chiral theory [31], minimally coupled 
fermions and gauged critical bosons, as well as minimally coupled bosons with gauged 
critical fermions 22, ^7], 26], see also [16] for a preliminary suggestion for duality. Our 
results may be viewed as additional evidence in support of these conjectured dualities. 

In order to find explicit formulas for the partition function of Chern-Simons matter 
theories on S 2 we need to evaluate the expectation value ( |1.3| ). This may be achieved using 
path integral techniques. Infact the partition function of Chern-Simons theories on S ff x S 1 
was evaluated by using path integral techniques long ago in the beautiful older paper |32| . 
the analysis |^| is easily generalized to include the effect of v(U). Briefly (see section [| 
for details) we follow [32] and work in the 'temporal' gauge d^A^ = 0, (in this paper the 
Euclideanized time direction is x 3 so ^3 denotes temporal component of the gauge field). 
We then abelianize the residual two dimensional gauge invariance by an additional gauge 
fixing condition; the holonomy matrix U(x) = e l ^ Aa ^ is required to be diagonal. The 
integral over the off diagonal and Kaluza Klein modes of the spatial part of the gauge field 
is quadratic and yields a determinant, which turns out to largely cancel the Fadeev-Popov 
determinant of gauge fixing. The integral of the diagonal (unfixed abelian part) of the two 
dimensional gauge field yields a delta function that fixes U{x) to be constant on S 2 . And 
the summation of U(1) N flux sectors discretizes the eigenvalues of U in units of ^x- These 
manipulations reduce the expectation value ( |1.2[ ) to a 'discretized' version of integral over 
the holonomy matrix U (|l.l|) 



N 00 



m=l n m =— 00 



n 2sin 



ai(n) - a m (n) \ -v(U) 



(1.6) 



where 



V(U) = T 2 V 2 v(U) 

and U(n) is the unitary matrix whose eigenvalues are e %am<yn \ where a m (n) 



(1.7) 

Z) with m running from 1 to N. Notice that the discretization interval between two allowed 



2irn„ 



values of eigenvalues in ( |l.6[) , tends to zero in the 't Hooft limit k —¥ 00, N 



00. 



It follows that the eigenvalue density functions p(a) in discretized one fll.q ) and non- 
discretized one (|1.1| ) obey identical large iV saddle point equations upon identifying the 
two potentials. Nonetheless the saddle points of ( |1.6|) and (1.1) in the 't Hooft limit are 
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not always identical. Recall that the classical large N saddle point descriptions of fll.6] ) 
and ( [l.lD are written in terms of the eigenvalue density function p(a) defined by 

1 N 

p( a ) = N^2 5 ( a ~ a ™)- ( L8 ) 

As the summation in ( |1.6| ) effectively excludes terms any two n m equal, 8 it follows that the 
maximum number of eigenvalues in an interval Aa in the summation in (|l,6| ) is given by 
4f, so that the eigenvalue density fll.8D is bounded from above by ^ x = ^j- In other 



words the eigenvalue distribution in (1J3) is constrained to obey the inequalities 



< P (a) < —. (1.9) 

On the other hand p(a) for the integral ( |1.1| ) obeys the lower bound listed in (|0]) (this 
is because a density is an intrinsically positive quantity) but no upper bound. 9 Note that 



the upper bound in (|1.9| ) agrees perfectly with the conclusions of |p7 |, obtained using 
Hamiltonian methods for Chern-Simons theory on T 2 . 

Any saddle point of (|l.l|) that happens to everywhere obey the inequality (|l.9|) is 
also saddle points of (|1.6| ). However saddle points of ( |Tl|) that anywhere violate the upper 



bound (|1.9| ) are not large N solutions of (|T^). Instead (|1.6|) admits new classes of solutions; 



those that saturate the upper bound of the inequality (1.9) over one or more arcs along 



the unit circle. While saddle points of (1^) may be classified in terms of the number of 



gaps in the solution (i.e. the number of arcs over which the eigenvalue density vanishes), 
saddle points of (|0]) are classified by the number of 'lower gaps' (regions over which the 
eigenvalue density vanishes) together with the number of 'upper gaps' (arcs over which the 



upper bound of (L9) are saturated). 

In order to understand the generic phase structure matrix model ( |1.6P we found it 
useful to first study a toy model in which V(U) takes the simple Gross -Witten -Wadia 
form V{U) = —-S(TrU + TiU^). We have found the exact solution to the saddle point 



equations of this toy model. The phase structure of this solution is depicted in Fig. |6(a) (see 



7\fy for a a more detailed summary of this solution); we pause to elaborate on this phase 



3 In the case that the base manifold is an S 2 , as studied in our paper, the exclusion of such configurations 



follows from the fact that the measure factor in (1.1) eliminates their contributions. The generalization 
of ([hi]) to the partition function of Chern-Simons theory on E 9 x S 1 where E 9 is a genus g manifold of 



arbitrary metric is given by the formula (4^6). In these cases the measure factors is either constant (in the 
case of g = 1) or diverges when two eigenvalues are equal. Nonetheless the correct prescription (the one 
that agrees with Chern-Simons computations using other techniques) appears to be to omit the contribution 
of such sectors. The justification for this prescription does not appear to be clearly understood from first 
principle path integral reasoning. We hope to clear up this point in the future. We thank O. Aharony, S. 
Giombi and J. Maldacena for extensive discussions on this point. 

9 The same upper bound for the density function appears in two dimensional Yang-Mills theory (p = 1), 
in which case the situation becomes similar due to the fact that there is no propagating degrees of freedom 
for the gauge field. A phase transition relevant to this upper bound was studied in 2d (q-deformed) Yang- 
Mills theory on S 2 |3l|, Ef^, which we will see from Chern-Simons theory below. We noticed these 
relevant papers when we were completing this paper. 
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diagram. In this toy model it turns out that for A < A c = ^ the eigenvalue distribution has 
no gaps for ( < 1, a single lower gap for 1 < ( < and a lower gap plus an upper gap for 
C > Try. This sequence of eigenvalue distributions is depicted graphically in Fig. 7(a) , 7(b) 



7(c), 7(d). For A > i, on the other hand, the eigenvalue distribution has no gaps for 



C < j — 1, a single upper gap in the range j — 1 < ( < 4X (i-\) , and an upper plus a 



lower gap for still larger £. This phase sequence is depicted in Figure 8(a) , 8(b) , 3(c) , |S(d) 



in §7.6. Moreover, in the large £ limit the eigenvalue distribution tends to the universal 



configuration 

p(a) = — !— (\a\ < 7rA) 

PV ; 2vrA yi 1 ; (l.io) 
= (|a|>7rA) 

at every value of A < 1. The distribution ( |1.10| ) is the nearest thing to a 5 function 
permitted by the effective Fermi statistics of the eigenvalues in perfect agreement with the 
results of the [27] which was the inspiration for the current work. 



For any given matter Chern-Simons theory V(U) is more complicated than the toy 
model of the previous paragraph. Even for these more complicated functions V(U), how- 
ever, it turns out to be easy to solve exactly for the eigenvalue density function in the 
'no gap' phase (see section || for details). The sequence of phases described above for the 
toy model is also true for every particular CS matter theory we have studied. Our exact 
solution allows us to determine the first phase transition temperature as a function of A as 
well as the value of A c in the theories we have studied (the value of A c varies from theory to 
theory). On general grounds we expect the eigenvalue distribution to tend to the universal 



function ( |1.1C ) in the high temperature limit in all Chern-Simons matter theories, even 
though we have not (yet) found exact solutions to the relevant matrix models in the two 
gap phase. 10 

We pause to describe an important property of our solution to the matrix model from 
real Chern-Simons theories in the low temperature no gap phase (for notational simplicity 
we work on a round sphere of volume Vi = 4ir in the rest of this paragraph). In the low 
temperature phase the partition function takes the form 



In every theory we have studied, the Taylor expansion of the function /(^) takes form 

f = a N2 +b{x) m + --- (L12) 



10 Note that the case A = 1 is special. In this case the number distinct allowed 'slots' for the eigenvalues 
(number of distinct allowed values of n m ) is exactly equal to N. As a consequence the eigenvalue distribution 
at A = 1 is given by 

P(«) = ^ (1-11) 

at all values of the temperature. In this strong coupling limit the holonomy distributes uniformly around 
the circle. 
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where a is a constant independent of A. It follows that at very low temperatures 

InZ = aT 4 (l + 0{T 4 /N 2 )) . 

This is precisely the form of the partition function of a gas of multi 'traces'; the value of the 
coefficient a also works out to the value predicted by this expectation. The fact that the 
free energy is independent of A at temperatures T 2 -C N is in perfect agreement with the 
non renormalization theorem of the spectrum of single trace operators in these theories. In 
other words the results for the partition function obtained in this paper at small -jr, match 
smoothly onto the high temperature expansion of the A independent O(N ) free energy. 

We have explained above that the eigenvalue distributions of two theories related by 
level rank type dualities are related by the map 

trU n (-l) n+1 tiU n . 



This relationship implies that the saddle point eigenvalue densities of two dual theories are 
related by the formula 

~ p{a) = 2^h\) ~ l^X p{a + n) 

where p(a) is the eigenvalue distribution of the theory with coupling A while p(a) is the 
eigenvalue distribution of the theory with coupling 1 — A. An interesting feature of this 
map is that it interchanges lower and upper gaps, consistent with the fact that the first 
phase transition is always from the no gap to the lower gap phase at small A while it is 
from the no gap to the upper gap phase at A near unity. 

We end this introduction on a cautionary note. All of the results of this paper have 
been obtained starting from the effective action (|1.2| ) . While the arguments for this effective 
action presented in [27] and reviewed in subsection ^T] below are persuasive, they are not 



completely beyond question. It would be useful to present a clear justification of (1.2) and 
understand how to compute corrections to the results from this action (such corrections are 
presumably suppressed in jt). It would also be useful to better understand the exclusion 
of equal eigenvalues from ( f4.6[) . We feel that it is important to clear up these two issues 
before we can regard the large number of exact results presented in this paper as beyond 
reproach. We hope to return to these issues in the future. 



2. Effective action and level rank duality 

In this section we review the effective description of large N Chern-Simons theories at 



temperature of order y/N proposed in [27]. This effective description allows us to relate 
the partition function of any fundamental matter Chern-Simons theory, in an appropriate 
high temperature limit, to the expectation value of an appropriate 'Wilson loop' in the 
pure Chern-Simons theory. Using this fact we then explore the implications of level rank 
duality of Wilson loops in pure Chern-Simons theories for the partitions under study. 
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2.1 The high temperature effective action 

Consider a Chem-Simons-fundamental matter theory on S 2 x S 1 in the 't Hooft large N 
limit. As described in the introduction, we study a theory at very high temperatures with 

V 2 T 2 = N( (2.1) 

where £ is held fixed in the large N limit. 

Our strategy for evaluating the partition function is, roughly speaking, first to find an 
effective action for the two dimensional holonomy field, U(x), around the thermal S 1 (here 
x is a point on the base S 2 ) by integrating out all other fields, and then to perform the 
path integral over two dimensional unitary matrix valued field U(x). In practice we find it 
more convenient to break up the first step (computation of the effective action of V(U)) 
into two steps as we now explain. 

The effective action for U(x) is obtained by summing all vacuum graphs in which all 
non- holonomy fields appear in an unrestricted manner. We will find it useful to break up 
these graphs into those that involve at least one matter field, and those that do not. Let 
the result of summing over all graphs of the first sort be denoted by e~ Se ff*- u ( x ". Once we 
have S e ff(U), the remaining task is to sum over all graphs that do not involve a matter 
field; and then integrate over U{x). 11 In other words the full partition function is given in 
terms of S e ff(U(x)) by the formula 

Z CS = J DAe^fi^+l^h^ffM. (2.3) 

While the formula ( |2.3[) is always correct in principle, it is useful in practice only when 



S e ff(U) is a local function. The key observation of [27] is that S e ff(U) is in fact ultralocal 
in the leading large N limit, when the temperature is scaled as in ( |2.1[ ). At leading order 
in large N, in other words the effective action 

S eff (U)=T 2 J d 2 x^v(U) (2.4) 

(the power of temperature has been inserted on dimensional grounds, see below for more 
details) and the formula for the partition function are simplified to 

Z cs = f DAe^K^+f^-^/^Vs »(C0. (2.5) 



The path integral in ( [2.5[ ) is evaluated in Euclidean space and the temporal circle is compact 
with circumference (5. 



11 Our conventions of the gauge field are those of i pTf . That is, the gauge field A M takes values of 
antihermitian matrix, the gauge covariant derivative is = 9 M -+- and its field strength is F M „ = 
[Dp,D v ] = d^A v — duAfj, + [Aft,A v ]. The gauge transformation is A^ — > A^ — [_D M ,e], where e is a gauge 
parameter taking values of antihermitian matrix. The Chern-Simons action is 

^ f Tr (MA + ^A 3 ^ =JLf d 3 xTre^ p (a^A p + ^A v a)j (2.2) 
where e 123 = 1. Under an infinitesimal change of A the action above changes as ^rTr J 5A A F. 
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The argument for the locality of <S e // is essentially that the graphs that we integrate 
out to find S e ff always include a matter propagator, and all matter fields develop large 
thermal masses. It appears to follow that S e ff(U) can be expanded in a series of local 
operators in the form 

S eff = J d 2 x (T 2 v(U) + Vl {U)TxDiUD l U + . . .) (2.6) 



[2 



where we have inserted powers of T by dimensional analysis. Now the scaling (2.1) 
converts the high temperature expansion (|2.6| ) into an expansion in inverse powers of N. 
For instance the terms listed above are respectively of order iV 2 , and N. 13 Terms that 



we have ignored in this expansion (the ... in (|2.6| )) are further suppressed at large N. At 



leading order in the large N expansion it seems to be justified to truncate (2.6) to the first 
term, yielding (|2.E|). 14 

The arguments leading up to ( |2.5| ) could conceivably have loopholes; in order to gain 
more confidence in ( |2.5D it would be useful to estimate the contributions of the neglected 
terms in (^1]) to the partition function and verify that their contribution is indeed sub 
leading in the jj expansion. We leave this exercise to future work. 

In order to make ( |2.5| ) a concrete formula we must explain how the function v(U) may 
be computed for any particular vector matter Chern Simons theory. We take this question 
up in the next section, and present exact and reasonably explicit formulas - as a function 
of the 't Hooft coupling A for v(U) for sample matter Chern-Simons theories. In the rest 



of this section, however, we first explore a structural 'duality' of the formula ( |2.5[) . 

2.2 Relations from level rank duality of pure Chern-Simons 

The equation (^) may be rewritten as 

Zc S = (e- T2 f**^ v W*») N , k (2.8) 

where the symbol 

denotes the expectation value of \I / in pure U(N) Chern-Simons theory at level k. As pure 
Chern-Simons theory is topological, an expectation value of the form 



/ e w(tf(s))' 



12 We are emphasizing that the S 2 partition function of fundamental matter -Chern-Simons theories are 
trivial (and A independent) at temperatures of order unity. These partition functions have interesting (and 
generically A dependent) structure only when the temperature is scaled as in ( |2.l[ ). 



13 Every term in the effective action (2.6) is of order N at fixed T (powers of N are counted as explained 
in the introduction), as the action is generated by integrating out fundamental fields. 

14 jijl considered a slightly different effective action, one in which they integrated out all graphs of matter 
fields, Kaluza Klein modes of the gauge fields, but not the massless 2d gauge field. They expanded their 
effective action in the form 



Opt 



J d 2 x (t 2 v w (U) + ui([/)TrA:E/D ,; l7 + v 2 (U)e ij TiFij + -j, « 3 (E/)TrF I3 F lJ + . . A ; (2.7) 

once again all terms but the first are sub leading in the j? expansion, ft would be interesting to carefully 
investigate the relationship between S e ff and Sj//- 
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is independent of x. It follows that (2.8) may be rewritten as 

Zcs = (e~ T2v ^) N , k . (2.9) 

Recall that v(U) is a gauge invariant function of the holonomy matrix. Any such function 
may be expanded in a basis of characters of U(N), i.e. in a basis of Polyakov or Wilson 
loops, in arbitrary representations of the gauge group U(N), around the curve that winds 
the thermal circle. In other words the remarkable formula (^1]) relates the thermal partition 
function of fundamental matter Chern-Simons theories (in an appropriate coordinated high 
temperature and large N limit) to the expectation value of a complicated thermal Wilson 
loop expectation value in pure Chern-Simons theory. 

Now pure Chern-Simons theories of rank N and renormalized level k are dual to pure 
Chern-Simons theories of rank k — N and level k via the well known level rank duality. 
Wilson loops transform under this duality as follows: a Wilson loop in the representation 
labeled by the Young Tableaux Y maps to the Wilson loop in the representation labeled by 
Young Tableaux Y where Y and Y are related by 'transposition' (the interchange of rows 
with columns). This duality map for Wilson loops may be restated in a slightly simpler 
fashion for arbitrary gauge invariant functions (like e v ^) of the holonomy matrix U as 
we now explain. 

The Wilson loop in any representation of U(N) is expressed as a polynomial of the 
variables TiU n (n < N) by the character polynomials of group theory. Let xy(U) denote 
the character polynomial of the U(N) representation with Young Tableaux Y. Let n denote 
the number of boxes in the Young Tableaux Y. The character polynomial of Xy(U) is given 
by the so called Schur Polynomial (see e.g. eq (2.15) of |38tl) 



= ^ E XrW ( ft (Tr£/ m ) fc "> j (2.10) 

a&Sn \m=l / 

where the summation on the RHS runs over all elements a of the permutation group S n , 
Xy(g) is the character of the permutation group element g in the representation labeled 
by the Young Tableaux Y and k\, k% ■ ■ -k n denote the number of cycles of length 1, 2, ... n 
in the conjugacy class of the permutation element g. We now recall that characters of the 
permutation group have a well known transformation property under transposition of the 
Young Tableaux 

Xy(9) = sgn(g)xY{g) (2.11) 

where sgn(g) is the sign of the permutation g (see e.g. [)39|j). The sign of the permutation 
g is given in terms of its number of cycles k n by 



-li k 

m=l 

Plugging (|2~TT1 ) and ( gig ) into ( plop we find 



sgn(g)=U [(-l) m+1 ] km . (2.12) 



X?M 4E Xy ^ ( A ((-ir +1 Trtn fc - ) • (2-13) 



n. 

ff£S„ \m=l 
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In other words transposition of Young Tableaux in all character polynomials (i.e. in arbi- 
trary gauge invariant functions of U) is achieved by the uniform interchange 15 



TrU n O (-l) n+1 Tr[/ n . 
It follows that level rank duality of pure Chern-Simons theory asserts that 



)k-N,k 



(2.15) 



(2.16) 



for any gauge invariant function V(U). 

In the large N limit, it is sometimes more convenient to regard V(U) as a functional 
of the eigenvalue distribution of p(a) than as a function of TrU n . Recall that TrU n = Np n 
where p n is the n th Fourier mode of the eigenvalue density function p(a). It follows that 
the interchange ( p. 15 ) amounts to 

N Pn = (k-N)(-l) n+1 p n . (2.17) 

Where p n are the Fourier modes of the holonomy in the CS theory with rank k. ( |2.17| ) 
implies 

\n+l N n _ /-gn+l ^_ 



Pn = (-iy 

(for n/0). For n = we define po 



~i \ \ Pn Ana 



k-N Pn " -' 1-A 
po = 1. It follows that 

1 A 



2tt(1 - A) 1-A 



p(a + 7r) 



or 



(1 - A)p(a) + \p{a + tt) = — . 
( 2.16 ) is thus equivalent to the assertion that 

V\o\\ 

)k-N,l~X 

where the functional V is defined via the relationship 

V[p] = V[p] 

for arbitrary functions p(a). 



(2.18) 

(2.19) 
(2.20) 

(2.21) 
(2.22) 



3 For the reader who dislikes appealing to mathematical authority for proofs, we present a 'physics' check 



of (2.13). Let Xk denote the character polynomial of the representation with k boxes in the first row of the 
Young Tableaux and no boxes in any other row. In a similar manner let tpk denote the character polynomial 
of the representation with k boxes in the first column of the Young Tableaux and no boxes in any other 
columns. Then it follows from the usual formulas of Bose and Fermi statistics that 



exp 



& fc2 

fe=0 
oo 

ipkX k = exp 



E 

m— 1 
oo 

E 



TrU m x r ' 



(-iy 



l TrU m x r ' 



in 



(2.14) 



where x is any real number. Level rank duality interchanges tpk and \k- At the level of the generating 
functions in (2.14) this interchange is simply achieved by (2.15). 
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2.3 Implications for Chern-Simons matter theories 

Let us label the class of fundamental matter Chern-Simons theories by an abstract label a. 
As we have emphasized above, the partition function of such theories at rank N and level 
k is given by an expression of the form 

where the index a emphasizes that depends on the particular Chern-Simons matter 

theory under study. 

Now suppose it happens to be true that for two different CS matter theories, theory a 
and theory b 

VZ x [p} = V'_ x [p}. (2.23) 



It then follows from ( |2.21 ) that the partition function of theory a at rank N and level k 



is identical to the partition function of theory b at rank k — N and level k. The identity 
of partition functions holds at all temperatures (of order y/N) on S 2 and also on arbitrary 
genus Riemann surfaces. 

In the next section we will find examples of simple Chern-Simons matter theories that 



obey (|2.23[) , strongly suggesting (previously conjectured) level rank type dualities between 
the relevant non-topological matter Chern-Simons theories. 

p£{ , |29[ ] had previously supplied evidence for several super symmetric Giveon-Kutasov 
theories by demonstrating that the S 3 partition functions of these theories - as computed 
by super symmetric localization - are in fact equal. In fact the equality of S 3 partition func- 



tions was demonstrated in [28, |25j by relating these partition functions to the expectation 
value of unknotted Wilson loops in pure Chern-Simons theory, and then demonstrating the 
equality of the resulting expressions using level rank duality of pure Chern-Simons theory. 
It is interesting that a similar mechanism appears to ensure the equality of thermal parti- 
tion functions of the dual theories (at least at high temperature), even though the thermal 
partition function (unlike the S 3 partition function) is far from a topological quantity. 
This suggests that it may be possible to prove duality between matter Chern-Simons the- 
ories starting from the level rank duality of pure Chern-Simons theory. We leave further 
development of this idea to future work. 

3. Computation of v (U) 

In the previous section we explained how the high temperature large N limit of the partition 
function for any fundamental matter Chern Simons theory takes the form ( |2.9[ ) for some 
function v(U). The form of the function v(U), of course, depends on the specific Chern- 
Simons fundamental matter theory under study. In this section we present exact results 
for the potential function v(U), as a function of the 't Hooft coupling A, for several simple 
Chern-Simons matter theories. We then go onto observe that these potential functions 



happen to be related by ( 2.23 ) for three pairs of matter Chern-Simons theories that had 



previously been conjectured to be related via level rank type dualities. Our results may 
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be regarded as further evidence for both these dualities as well as the correctness of the 
effective description (JO)|), 

We have obtained our results for v(U) by making only minor modifications to several 
existing computations in the literature. In order to explain this we pause to present a 
brief review of the recent relevant literature on fundamental matter Chern-Simons theories. 
Such theories at nonzero A have recently received serious attention, starting with the papers 



1 16|, [L8|] and subsequently [^T], 23 , |4|] • I n particular the authors of |L6| worked in a lightcone 
gauge to present an exact, all orders computation of the infinite volume finite temperature 
partition function of the fundamental fermion Chern-Simons theory in the background of 
the identity holonomy matrix. 16 The recent paper [27] pointed out that the computation of 
fl6fj is easily generalized to an other holonomy background; |27| repeated the computation of 



PI (and generalizations) in the holonomy background ( l.lOj ). In fact the partition function 



computation of [16] may easily be repeated in an arbitrary holonomy background. In this 
paper we propose that the result of this computation may be identified with T 2 V 2 v(U). 
This proposal allows us to explicitly compute v (U). The computations that yield v(U) are 
a straightforward generalization of the work of [23|, 



| 27 1 . For this reason we do not 
describe the method employed to evaluate v(U) but simply present our final results below. 



3.1 CS theory coupled to fundamental bosons with classically marginal inter- 
actions 

In this subsection we study the Chern-Simons matter theory coupled to massless funda- 
mental bosons with a 6 interaction. The Lagrangian of the theory we study is presented 
in equation (4.2) of |27]. 
We find 

at o /V f n f°° 

v[p] = -TT^i 1 + + o - / da P^ / d V y N 1 " e ~" +iQ ) + ln ^ " e ~ V ^) < 

6vr A(A,A 6 ) 2vr J a 

(3.1) 



where 



A(A,A 6 



(3.2) 



The potential v[p] above is determined in terms of the constant a that has a simple inter- 
pretation; it determines the thermal mass of the bosonic field via the equation 



a 2 T 2 . 



(3.3) 



16 The authors of Jl(| (and extensions to other theories [jL6|, p3|, |3(]]) proposed that the partition function 
evaluated in these works should be identified with the thermal partition function of the relevant theories. 
This identification was later realized to be inconsistent |23),j3^] with Giveon-Kutasov duality for a chiral 
theory The very beautiful recent paper |2^] pointed out that this identification p3l ^| is incorrect. 
They used the Hamiltonian methods of |}40(] to explain that the R 2 thermal partition function of matter 



Chern-Simons theories is given by the partition function of the relevant theories in the holonomy (1.10) 



This proposal yields results consistent with Giveon Kutasov duality, as well as a new non supersymmetric 
version of this duality conjectured and tested in Ji^, 
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The value of a is determined by the requirement that it extremizes v[p] at constant p, i.e. 
a is required to obey the equation 



1 A^6 7^, r , i \ ( ■> r. • , ,o — ioi . .a + ia. 
a = — \ — \ + \ 2 l dap(a) I In 2 smh( ) + ln2smh( 

2 V 8ir J— n V * 



2 ' (3-4) 
b> - J da\p(a) (in 2 sinh( — - — ) + In 2 sinh( — - — ) J 



A 



The expressions presented above agree with previous computations in the literature for 
special choices of p(a). Note in particular that, (|3.4|) reduces to Eq.(4.22) of (27j upon 
replacing the arbitrary density function in (3.4) with fll.10 ). Similarly (3.1) reduces to 
Eq.(4.33) of p7f] upon making the same replacement. 

3.2 CS theory minimally coupled to Fermions 

In this subsection we study the Chern-Simons matter theory coupled to massless funda- 
mental fermions. The Lagrangian of the theory we study is presented in equation (2.1) of 
We find 

v[p] = -^-^j-c 3 + 3j 7r dap(a) J™ dy y(ln(l + e^"*) + ln(l + e~^))) , (3.5) 

where c once again determines the thermal mass of the fermions. 17 The value of c is 
obtained by extremizing v[p] w.r.t c at constant p, i.e. c obeys the equation 



A 



J* da p(a) ^ln 2 C osh(^^) + In 2 cosh(^^) J . (3.7) 



The expressions presented above agree with previous computations in the literature for 
special choices of p(a). Upon setting p{a) = 5(a), (^) and ( |3.7|) reduce to Eq.(2.80), 
Eq.(2.74) of || while fiufy reduces to Eq.(2.92) of p|. On the other hand, on setting 
p(a) to the expression given in ( 1.1CQ , ( [3~7| ) and ( |3.6|) reduces to Eq.(5.13), Eq.(5.14) and 
Eq.(5.15) of |^7| (upon setting a = where a is defined in that paper ) while (|3.5|) reduces 
to Eq.(5.28) of [E7| in appropriate limit as discussed above. 



3.3 Chern-Simons coupled critical boson 

In this subsection we study the Chern-Simons matter theory coupled to massless criti- 
cal bosons in the fundamental representation, i.e. a Chern-Simons gauged version of the 



17 As in the bosonic case Er = c 2 T 2 is the thermal mass of the fundamental fermions. More precisely, 
the fermionic self energy is given by £r(p) = f(/3p a )p 3 I + ip~ g(/3p 3 ) r y~ , where 



f(y) = - / dap(a) ln2cosh(^ ) +m 2cosh(^ ) 

VU V 2 2 / (3.6) 

C 2 2 
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U (N) Wilson Fisher theory. The Lagrangian of the UV theory is simply that of massless 
minimally coupled fundamental bosons deformed by the interaction 



5S 



(3.8) 



where A is a Lagrange multiplier field (see subsection 4.3 of p7| for details, and in particular 
eq.(4.35) for the Lagrangian; note |27]] employs the symbol a for our field A). We find 



v[p\ 



ydy I p(a)da (ln(l - e ~ y+ia ) + ln(l - e~ y - ia )) 

(3.9) 



The squared thermal mass of our system is a 2 T 2 . The constants A and a above are 
determined by the requirement that they extremize v[p] at constant p. This requirement 
yields the equations 

(3.10) 



A = a 2 T 2 , 



and 

r ( ■ a 

2 ' v 2 

The expression for v[p] simplifies somewhat upon inserting ( 3.1C| ) into fl3,9|) ; we find 



p(a) ln2sinh( )+ln2sinh( — - — )) = 0. 



(3.11) 



v[p] 



-^ 3 + ^ I >'!' I .W(n)(ln(l )+ln(] -.-"-'")). (3.12) 



where a is obtained by extremizing ( 3.12j ), i.e. it is determined by the equation ( 3.1 1| ) . 
Notice that neither ( |3.12 ) nor ( p.ll ) depend on the 't Hooft coupling A. It follows that v[p] 
is independent of A as well as A6 (this was not the case in the previous two subsections). 

The expressions presented above reduce to formulas that have previously appeared in 
the literature for special values of p. In particular for p given by ( l.lOj ), ( |3.12 ) reduces to 
Eq.(4.41) of |§, while ( |3lID reduces to Eq.(4.42) of (2^]. 



3.4 Chern-Simons theory coupled to critical fermions 

In this subsection we study the Chern-Simons matter theory coupled to massless critical 
fermions in the fundamental representation. The Lagrangian of our system is simply that 
of massless fundamental fermions minimally coupled to Chern-Simons theory deformed by 
the terms 

(3.13) 



5C = BinP + ~—\ f 6 B 3 



N 
~6 



where B is a Lagrange multiplier field and Ag is a coupling that is marginal in the critical 
fermion theory as B turns out to have unit dimension. We find 



v[p] 



N f°° , 

SF/, * 



l 



3^ + (^_ , 



2A 



2A 



(3.14) 



da y p(a)(ln(l + e 



y-ia\ 



+ ln(l + e 



-y+ia 



)), 



-17- 



where c T is the thermal mass of the fermions. 18 The constants B and c above are deter- 
mined by the requirement that they extremize v[p\. This requirement yields the equations 



(3B = ± 



aa /3(a) ln2cosh( — - — ) + m2cosh( — - — 



(3.16) 



1 - 27TAA, 



Assuming that c is positive, we are forced to choose the negative root in the second of 
( 3. 16 ). 19 Inserting the second equation in ( 3.16 ) into the first we find 



c 1 + 



1 



A 



1 - 2vrAA 



/ 



da p(a) ( ln2cosh( C ^/ a ) +ln2cosh( C %a ) ) . (3.18) 



Inserting the same equation into (3.14) yields 
N 



v[p\ 



where 



6vrA 



1-A + $(A,A£))+3A 



,hj j p(a)da y(ln(l + e~ y - ia ) + ln(l + e~ y+ia )) ) , 

(3.19) 

1 



1 - 2vrAA 



(3.20) 



The expressions presented above reduce to formulas that have previously appeared in the 
literature for special values of p. In particular for p given by ([LT0|) , ( [3l8l) reduces to 
Eq.(5.35) of p7|] , while ( |3.19|) reduces to Eq.(5.34) of 27] with appropriate factors. 

3.5 M = 2 theory with a single fundamental chiral multiplet 

The Lagrangian of the theory we study is presented in equation (2.1), (3.58) of [^] (see 
also (6.1) and (6.20) of pi). We find 



N 



v[p\ 



6vr|A| 



c 3 — 6 1 A | / dap(a)Ke I (iyylogtanh 



y + ia 



(3.21) 



18 The fermion self energy (the correction to the quadratic term in the fermion effective action) is given 
by S T (p) +BI — f(Pp s )psI + ip~ g{l3ps)-y~ , where 



X 



f(y) = -B + - / da p(a) ln2cosh( 
y V 



y 



-y/j/ 2 + C 2 + ! 



+ In 2 cosh ( 



(3.15) 



c — /3B = A / da p(a) ( In 2 cosh(^-^) + In 2 cosh(^— — ) 



19 This may be seen as follows. Inserting the second equation in (3.16) into the first equation of (3.16) 
we obtain 



c it 



1 - 2ttAA 



/ 



A / da p(a) I In 2 cosh( ° — iQ ) + In 2 cosh(- — — 



(3.17) 



As the R.H.S. of this equation is positive, we must choose the negative root, i.e. j3B = 
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where for the thermal mass (for both the boson as well as the fermion in the supermultiplet) 
is denoted by cT. The constant c above is determined by the requirement that it extremizes 
v[p}\ i.e. 

c + ia 



/7T 
da\p(a) log coth ■ 
-7T 



- (3-22) 

The expressions presented above reduce to formulas that have previously appeared in the 
literature for special values of p. In particular for p given by ( 1.10Q , ( |3.22[ ) reduces to 



Eq.(6.22) of pfl] , while (|3.2l| ) reduces to Eq.(6.23) of [27] with appropriate factors. For the 



case when p(a) = 6(a), (|3.22|) , ( 3.2l| ) are straight-forward generalization of (3.61), (3.62) 

of m. 



3.6 Level-Rank duality 

The Chern-Simons matter theories we have studied, earlier in this section, have been 



conjectured to be related to one another by strong weak coupling dualities [20, |2^, [27L 26], 



see also [16] for a preliminary suggestion. The best established of these conjectures is the 
conjectured strong weak coupling self duality of the supersymmetric theory. According 
to this conjecture the supersymmetric theory at rank N and renormalized level k may 
be identified with the same theory at rank k — N and level k. The non supersymmetric 
theories we have studied above have also been conjectured to be related to one another by 
strong weak coupling type dualities. In particular the theory of minimally coupled fermions 
at rank iV and renormalized level k has been conjectured to be identical to the theory of 
Chern-Simons gauged critical bosons at rank k — N and level k. 

Finally, the minimally coupled scalar theory and the critical fermion theory are con- 
jectured to obey the following duality. The minimally coupled scalar with rank N, 't Hooft 
coupling A and A(A, Xq) is conjectured to be dual to the fermionic theory with rank k — N, 
t' Hooft coupling 1 — A and g(l — A,Ag) = - ^ - (the last relationship effectively deter- 
mined the relationship between the two marginal parameters) where A, g are defined in 
(|3^),(^20|) respectively. 

In this subsection we will demonstrate that the potentials, v[p], for (conjectured) dual 
pairs all obey ( |2.23 ). It follows that all the partition functions computed in this paper 



are consistent with all three conjectured dualities listed above. This agreement may be 
regarded as evidence in favor of all three dualities, as well as the effective action framework 
of section |2| on which the current paper is based. 

3.6.1 Self-duality of M = 2 



v[p] for the supersymmetric theory was reported in subsection 3J3. This potential depends 
on the values of N and A; in the rest of this subsubsection we use notation for this potential 
that makes its dependence on parameters explicit; we denote the potential by v s ^ s l[p\. We 
also denote the thermal mass of the supersymmetric theory by c™ sy [/o]. In the rest of this 
subsubsection we will demonstrate that 

<?rip] = ™M (3 23) 

susy r l susy r ~i V ' / 

V NllP\= V k-NA-\[P\ 
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where p is the transformed eigenvalue distribution function defined by the equation 

1 



(1 - X)p(a) = 2^- Xp{a + vr). 
Let us first demonstrate the first equation in (|3.23|) . According to ( |3.22| ) 

-susy r 



cTJ[p} = 2 



Re / da{l — X)p{a) log coth 



-l-AL 



+ ia 



Using (|3.24[) the RHS of fl3~22l) may be rewritten as 



Re / da ( Xp(a + tt) ] log coth 



+ ia 






= 2 


Re / 






Jo 



"2tt 



Re / daXp(a) log coth 



(3.24) 



(3.25) 



2 

(3.26) 



where we have used the integral identities 

x + ia 



da log coth ■ 



0, coth( — - — ) = tanh( — ) for all x. 



(3.27) 



It follows that 



c 



l-AL 



p2n 

Re / daXp(a) log coth 
Jo 



Re / da\p{a) log coth 



+ ia 



(3.28) 



Comparing with (|3.22j ) we see that c^^[p] obeys the same equation as c^ usy [/?]. Hence we 
conclude 



'l-AL 



demonstrating the first equation in ( 3.23 ). 

We now establish the second equation in ( |3.23 ). According to ( [3.21 ) 



susy 
' J k-N 



1-xlP] 



k-N 



(c?-l[p\) 3 ~ 6|1 - A| / dap(a)Re 



6 



6 



6 



6 



6 



A|vr|l - A| 
N 



y + ia 



dyy log tanh 

^-A Id 2 



A|vr 
N 



(5 



6 / da 



(CW-6 / da(-\X\p(a))Re / 

JO Jc- 



2tt 



\X\p(a + 7r) Re 



01 



(iyy log tanh 



y + ia 



X\tt 
N 
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(3.30) 



where we used ( p~27[) in going from the second to the third line, and ( |3.29| ) in going from 
the fourth to the fifth line. This completes our demonstration of the self duality of v[p] for 
the super symmetric theory. 
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3.6.2 Critical boson v.s. Regular fermion 



The potentials v[p] for the critical boson and the minimally coupled fermion theory were 



reported in subsections 3^3 and 3J2 respectively. v[p] depends on the values of N, A; in the 
rest of this subsubsection we use notation for this potential that makes its dependence on 
parameters explicit; we denote the potential for the critical boson theory and the regular 
fermion theory respectively by f^^'f/o], u ]vA F [p\- We also denote the thermal mass of the 
critical bosonic theory and the regular fermionic theory by <7^ n ' B '[/9], c^" eg,F '[/9] respectively. 
In the rest of this subsubsection we will demonstrate that 



v k-N,l-\[ 
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[p], 
[p], 



(3.31) 



where p is the transformed eigenvalue distribution function defined by the equation ( 3.24| ). 
To start with, we first demonstrate the second equation in ( 3.31[ ). According to ( 3.1 1| ) , 
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(3.32) 



Rearranging terms in this equation we find 
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Using the integral identity 
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( 3.33 ) may be rewritten as 
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Comparing ( 3.35| ) with 
establishing that 



we conclude that [p] obeys the same equation as c^ eg ' r ' [p] , 
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Now we proceed to demonstrate first equation of ( 3.11 ) . According to ( 3.12; ) 
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(3.37) 



where we have used ( |3.36| ) in going from the first to the second line. Comparing with (| 
we conclude 

v^ii-x[p]=v^f-[p}. (3.38) 
This completes our demonstration of ( p. 31 ). 



3.6.3 Regular boson v.s. Critical fermion 

The potential v[p] for the marginally boson and the critically coupled fermion theory were 



reported in subsections 3.1 and pA respectively. This potential depends on the values 



of JV, A, A6 for the boson, and on N, A and Ag for the fermion. In the rest of this 
subsubsection we use notation for this potential that makes its dependence on parameters 
explicit; we denote the potential of the regular bosonic theory and critical fermionic theory 
by v?ff*x '[p], v Gri ' F ' f [p] respectively. We also denote the thermal mass of the regular 

bosonic theory and critical fermionic theory by o^ e ^ ' '[p], c c "y F '[p] respectively. In the rest 
of this subsubsection we will demonstrate that 



6 



Consequently the relations ( 3.39| ) assert the duality between bosonic and fermionic theories 



with marginal deformation parameters related in a known but complicated manner. In the 
equation above p is the transformed eigenvalue distribution function defined by the equation 



d3~24| ) and is given as a function of Xq by the complicated relationship 
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2a/1-2ttAA^. (3.40) 
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Let us first demonstrate the second equation in fl3.39|) . According to ( fO| ) 
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(3.41) 



In going from the third to the fourth line of ( |3.4l| ) we have used the integral identity ( 3.34| ) . 
Rearranging terms, we obtain 
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Using (|3.40[) we obtain 20 
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Comparing ( 3.43j ) with ( 3.18 ) we conclude 
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fl3.40| ) is just a slight rewriting of Eq.(5.37) of Q which reads ^\f^z + Ag = 2^1 - 2tt(1 - A 6 )A£. 
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The first of the ( jOsD can be proved as follows. Under the duality one obtains 

Reg.B. r-i 
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Now using ( |QoD and the fact that da (ln(l - e"^ +ia ) + ln(l - e~ y -' ia )) = 0, one ob- 
tains 
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Using ( |3.20| ),( |3~T8| ) and ( 3.44 ) we conclude that above equation is same as ( |3.19| ). So we 
have shown 

^-A,A 6 [P]=^£/[P]- (3-47) 



4. Evaluation of the path integral ( |2.5| ) 



At first sight ( |2.5| ) is the path integral over a nontrivial interacting field theory in three 
dimensions, and its evaluation appears to be a daunting task. However the rewritten form 



( |2.8| ) emphasizes that (2.5) evaluates a topological observable in the topological Chern- 
Simons theory. This special features of Chern-Simons theory make the evaluation of the 
path integral almost trivial. The path integral we need was in fact largely evaluated in p2[ . 
In order to make our paper self contained, however, we reevaluate it below, very closely 
following |52]. Our focus in this section, and the rest of this paper, will largely be on 
Chern-Simons theories on S 2 x S 1 . However we will also briefly consider the generalization 
to Chern-Simons theories on T, g x S , where £ 5 is a genus g Riemann surface. 
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4.1 Gauge fixing 

We impose the following gauge fixing conditions: 

(i) d 3 U = dse^ 3 = 0. 

(ii) We simultaneously diagonalize U(x) for all points x on S 2 . 

(iii) We impose the Coulomb gauge on the time independent diagonal elements of Ai,A2- 

Here A3 is the gauge field along the S 1 direction and Ai^ are the ones along the S 2 
directions. 

The choice (i) is the nearest one that can come to the temporal gauge A3 = when 
the thermal time is compact. The gauge condition (i) leaves two dimensional gauge trans- 
formations (those that depend only on the coordinates of the sphere) unfixed. The gauge 
condition (ii) fixes some of this freedom, leaving only a two dimensional U(1) N gauge sym- 
metry unfixed. The gauge condition (iii) fixes the residual 2d Abelian gauge invariance by 
imposing (Abelian) Coulomb gauge for the time independent parts of the diagonal elements 
of Ai, where = 1, 2 are the sphere spatial indices. 

The gauge condition (ii) that we have employed to abelianize the residual two dimen- 
sional gauge invariance is very similar to the gauge employed by 't Hooft in his classic 
study of confinement due to monopole condensation. As noted by 't Hooft the gauge 
condition (ii) is ambiguous (i.e. fails to completely fix gauge) when two or more eigenvalues 
of the matrix U = e^ As coincide. Two eigenvalues of a unitary matrix coincide on a surface 
of codimension 3 in the space of all unitary matrices. In 't Hooft's study of confinement of 4 
dimensional gauge theories, such coincidences typically occurred on a 4 — 3 = 1 dimensional 
surface in spacetime, i.e. on a line. The world line of such coincidences were monopoles for 
the gauge fixed abelian theory; 't Hooft argued that the condensation of such monopoles 
is responsible for confinement. 

In the current paper the 'abelianizing' gauge condition (ii) is adopted on a two di- 
mensional manifold (i.e. the base manifold S 2 ). As a consequence we would generically 
expect eigenvalues to coincide ona2 — 3 = — 1 dimensional surface of our two dimensional 
base manifold (S 2 ); in other words in the two dimensional context of this paper two eigen- 
values of U(x) do not generically coincide. Consequently configurations with coincident 
eigenvalues do not represent 'objects' or 'defects' in our path integral. Nonetheless such 
configurations play an important role in our path integral as we now explain. 

Over the next two paragraphs we characterize the space of smooth functions U{x) 
before we impose gauge condition (ii). Smooth connections on S 2 x S l are all topologically 
connected to A^ = 21 ; in particular the set of all smooth functions U(x) form a single 
connected space. However our particular choice of gauge (ii) is singular when two eigen- 
values of the matrix U coincide. Such a singularity occurs on a subspace of codimension 1 
in the space of all smooth functions U(x) 22 . These singular configurations (as far as our 

21 U(N) connections have distinct topological sectors, labeled by the overall U(l) flux. 
22 This may be understood as follows. Consider a one parameter set of functions U(c,x) where c is a real 
valued parameter. Requiring that two of the eigenvalues of U(c,x) coincide give a set of three equations, 
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choice of gauge is concerned) constitute 'domain walls' that chop up the space of smooth 
SU(N) connections into a set of domains or unit cells. By definition, within a unit cell in 
the space of smooth functions U(x), the eigenvalues of U coincide nowhere on S 2 . How- 
ever on the domain walls two eigenvalues of U coincide somewhere on S 2 (more than two 
eigenvalues coincide at a point - or two eigenvalues coincide at multiple points on S 2 where 
domain walls intersect). As our gauge fails when this happens, we are forced to perform 
our path integral in each unit cell separately. Provided the path integral over every unit 
cell has no divergence when we approach the domain wall surfaces we can simply ignore 
the contribution of the domain walls to the path integral (as they are of codimension unity 
in field space) ; the full path integral is given by the sum of the path integrals in every unit 
cell. 23 

In order to see how this works in more detail, we now focus on the special case N = 2, 
the generalization to the SU(N) theory is straightforward. Within any unit cell U(x) is 
nowhere equal to ±7 where / is the identity matrix (this is true by definition of a unit cell). 
Within any unit cell, therefore, the function U(x) is a map from S 2 to S 3 where S 3 is the 
SU{2) group manifold S 3 with two points U = ±7 removed (the two points may be thought 
of as the 'north pole' and 'south pole' of the S 3 ). While maps from S 2 to S 3 are trivial, 
maps from S 2 to S 3 are characterized by an integer valued winding number. 24 We have 
reached an important conclusion; unit cells in the space of smooth connections are labeled 
by an integer valued winding number. Distinct unit cells have distinct winding numbers. 
The winding number is ill defined on the domain walls that constitute the boundaries of 
unit cells; i.e. at the points at which the gauge fixing (ii) will fail. 

We now turn to the implementation of the gauge condition (ii) (we continue to work 
with the special case of SU(2)). We wish to construct gauge transformation g(x) that 
implements the gauge fixing condition (ii); we proceed as follows. We first perform a 
global gauge transformation that diagonalizes U at the south pole of the base manifold S 2 
(note that the north and south pole of the spatial S 2 are completely distinct from the north 
and south pole of the group manifold S 3 discussed in the previous paragraph). We then 
diagonalize U(x) everywhere else on the S 2 as follows. Starting from the north pole, we 
move along lines of constant longitude (i.e. constant (j) using the usual conventions for polar 
coordinates on S 2 ) and determine a g(9,(fi) that diagonalizes U(x) along that longitude; it 
is easy to check that this is always possible. In the neighborhood of the south pole this 
construction gives us a function H{4>) = lim e _>o — e, <^>) , i.e. a map from a circle to the 

that will generically yield discrete solutions for the 3 variables c, xi, xi. As this condition is satisfied only 
at discrete values of the parameter c, it follows that the subspace of U (x) with two equal eigenvalues is of 
codimension unity in the space of all smooth functions U(x). 
23 We thank Xi Yin for a very useful discussion on this topic. 

24 This is analogous to the easily visualized fact that maps from S 1 to S 2 are characterized by an integer 

valued winding number, where S 2 is the sphere with north and south pole removed. Note also that for 

configurations for which A3 is small the winding number may be characterized in a slightly more elementary 

manner. Consider a smooth configuration Az{x) before imposing the gauge choice (ii). As long as As(x) is 

everywhere nonzero (i.e. provided its eigenvalues nowhere coincide), the unit vector field A%(x) = 7==g 

V 3 ^3 

is well defined on S 2 . This unit vector field is characterized by an integer valued topological invariant; it's 
wrapping number over the S 2 . 
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gauge group. Recall, however, that U was already diagonal at the south pole. In order 
that U continue to be diagonal after the gauge transformation it must be that H(<p) lies 
in the U(l) subgroup of SU{2) that leaves any diagonal element invariant. In other words 
H(4>) defines a map from S* 1 to U(l); such maps are characterized by an integer valued 
winding number. It is a famous result in the study of non abelian solitons that the winding 
number of the map H(<f>) is the same as the wrapping number of the map from S 2 to S* 3 
defined by U(x). 25 

Clearly the function H((f>) can be chosen to be a constant (so that the gauge transfor- 
mation g(6, is well defined on the south pole) only if its winding number vanishes. It 
follows that the gauge choice (ii) can be implemented by a smooth gauge function g(x) if 
and only if the wrapping number of U (x) on S 2 vanishes. 

When the wrapping number of U(x) is non vanishing on the S 2 , the singularity of the 
gauge transformation function at the south pole has a straightforward physical interpreta- 
tion. The winding number of the function H(4>) is simply equal, by Stokes theorem, to the 
U(l) flux that runs through the south pole. In other words once the gauge condition (ii) 
has been implemented starting with a function U(x) that had wrapping number n on the 
S 2 , the background U(l) connection has n units of flux over the S 2 (and a compensating 
Dirac string threading the south pole). Consequently when we implement our gauge con- 
dition (ii) within a unit cell characterized by winding number n, we are finally left with the 
path integral over the abelian theory on a sphere with n units of U(l) flux for the abelian 
gauge group. 

Let us summarize. Before fixing the gauge condition (ii) we were required to perform 
a path integral over each unit cell, labeled by a wrapping number of U(x) on S 2 . Upon 
fixing the gauge condition (ii) sum of the path integrals over each unit cell turns into a 
sum over path integrals in the abelian gauge fixed theory on an S 2 , where the abelian path 
integral is performed over connections with a nonzero abelian flux on the S 2 , equal to the 
wrapping number of U(x) on S 2 before we implemented the gauge condition (ii). In each 
flux sector we must now integrate unrestrictedly over the diagonal elements of U(x) and 
the abelian gauge fields A fl (x). 

In our presentation it might appear that the path integral in nonzero flux sectors has 
to be performed about a singular background, but of course the singularity of the gauge 
field at the south pole of S 2 is of the Dirac string type and is a gauge artifact. Within 
the residual U(l) theory we can desingularize this configuration a la Wu and Yang using 
coordinate patches. The path integral to be performed is perfectly smooth in every sector; 
it merely has to be performed in the background of smooth fluxes for the residual abelian 
group. We will implement this procedure in what follows. 



25 The function g{6,<j>) is not unique; if g(8,<f>) diagonalizes U(x) then the same is true of h(8,4>)g(8,4>) 
where h is an arbitrary U(l) diagonal gauge transformation. This ambiguity maps H((j>) to h(n,<f>)H (cf>). 
However as h(0, <j>) is independent of <f> (as g is chosen to be well defined at the north pole) it follows that 
h(-K, 4>) is smoothly connected to a constant (by varying 8) and so has no winding number. It follows that 
the winding number of H(<f>) is unaffected by this ambiguity and is also the only convention independent 
information in the function H(<f>). 
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4.2 Fadeev-Popov Determinants 

The Fadeev-Popov determinants for the gauge fixing conditions (i) and (ii) is 



A FP (A 3 ) = Det' s (D 3 ) = Det' s (d 3 + [A 3 ,]) (4.1) 

where Det' s is the determinant of the operator acting on all (monopole) scalar functions 
on the S 2 , in the given flux sector, exempting diagonal time independent modes. The 
gauge fixing condition (iii) fixes the residual abelian gauge invariance. As a consequence 
the associated Fadeev-Popov determinant is (divergent) constant, which will be regularized 
in a suitable way. In summary the full Fadeev-Popov determinant for our problem is given 
by (4.1). In Appendix [A| we heuristically evaluate and interpret this determinant. 



4.3 Evaluation of the path integral over vectors 

In our gauge the path integral fl2.5| ) reduces to 



Z = J DA 3 A FP (A 3 ) exp (-T 2 j d 2 x^g v(U)\ Z'[A 3 ] (4.2) 

where 

Z'[A 3 ] = J DA ie -^ +s ^ 
k? I' 

51 = — J d 3 xTr {A 3 (d x A 2 - d 2 A 1 )) (4.3) 

ki I' 

5 2 = — y d 3 xTv(D 3 A 1 A 2 ). 

It is easily verified that the action S± receives contributions only from time independent 
diagonal modes of Ai. It may also be verified that these modes (time independent diagonal 
modes of Ai) do not contribute to S 2 . As a consequence the path integral in ( [4.3| ) splits into 
the product of two path integrals, the first over time independent diagonal modes of Ai, 
and the second over all other modes. The second path integral is quadratic; its evaluation 
yields the determinant 

1 



^/Det' v (d 3 + [A 3 ,}) 

where by Det' v we mean the determinant of the operator D 3 acting in wedge product 
over (monopole) vector functions (i.e. the determinant of the operator in the quadratic 
form f (fixe 1 ! AiD 3 Aj) on the space of all (monopole) vector fields baring time independent 
diagonal modes. Note that this determinant generically depends on the flux sector we work 
in. 

We use the notation i am ^ (m = 1, . . . , N) for the N eigenvalues of ^3 and a l m {x) for 
the time independent diagonal modes of Aj. We work in the flux sector in which we have 
M m units of flux for the m th U(l) factor. The path integral over a l m (x) may be evaluated 
as follows. In each U(l) m factor, the most general gauge field configuration is given by a 
constant flux whose integral is given by 2irM m plus a fluctuating component of a l m that is 
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everywhere well defined on the sphere. The most general well defined vector field on the 
sphere is given in terms of two scalars by the formula 

(4(x) = SVmOc) + j j djXm(x). (4.4) 

The Coulomb gauge condition sets ip m (x) to zero. Substituting this solution into 52 in 
Q4.3D yields a term proportional to 



where V 2 is the Laplacian defined on the two sphere. The integral over Xm yields the delta 
function 

N 

11 -5(V 2 a m ). 

m=l 

This 5 function reduces the functional integral over a m (x) to the ordinary integral over 
the constant pieces a m . Consequently, the path integral over diagonal time independent 
modes in the relevant flux sector is simply given by 

da m e ikMmam . 

In order to evaluate the full path integral we must sum over all flux sectors. Before we can 
do that, however, we need to evaluate the dependence of the ratio of determinants on the 
flux sectors. We turn to that evaluation now. 

4.4 Evaluation of the ratio of determinants 

The ratio of the ghost determinant (our Fadeev-Popov determinant) and the determinant 
obtained from integrating over vector fields was very carefully performed in the paper of 
1 32], and we will simply quote their result below. In order to give the reader a flavor of the 
result, however, we present a hand waving 'evaluation' of the ratio of determinant (2.5) in 
the zero flux sector. 

Let us consider the vector quadratic form in more detail. Let us first work in the zero 
flux sector. We have 

' d 2 xe l3 A i D 3 A j . 

In the zero flux sector an arbitrary vector field can be written as di4> + eijd 3 x- Plugging 
in, we find that the quadratic form above is proportional to 

d 2 x (0Z? 3 V 2 x) 

where we have integrated by parts and used the fact that A% is a constant matrix. Ignor- 
ing possible subtleties involving regulation of these infinite products, it follows that the 
quadratic path integral over vectors is given by 

1 1 



y/Det' v {D z ) Det'^(D 3 
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where Det" s (D%) is the determinant of the operator D% acting on scalar functions ( |4.1| ) 
excluding the zero mode of the scalar field on S 2 . It follows that the ratio of the ghost 
determinant fl4.1|) and the determinant that arises from integrating over vector fields is 

Pet's ( D 3) _ TT 9 ■ ( am(n m ) - ai{ni) \ 



where the nonzero contribution to ( f4.5|) arises purely from the zero modes on S 2 (see 
Appendix |A| for the explicit formula). 



The result (4.5) was obtained much more carefully in the paper [32]. In that paper 
it was also demonstrated that ( |4.5| ) is valid on a space with the topology of the sphere 
independent of both the background metric of the space as well as the background fluxes 
(the Chern classes). We will not attempt to re-derive this result here, but will use it in 
what follows. 26 

The evaluation of the ratio of determinants presented here was heuristic because we 
cancelled infinite ratios of determinants without regulating them. It turns out that our 
heuristic derivation is actually accurate in the dimensional reduction regulation scheme 
which we use throughout this paper. In other regulation schemes, like regulation with a 
Yang Mills term - the determinant ratio has a phase which shifts the Chern-Simons level. 
Throughout this paper k refers to the renormalized Chern-Simons level; the effective level 
of the theory after accounting for this shift. 

4.5 Summation over flux sectors 

As the ratio of determinants is independent of the flux sector, the summation over fluxes 
operates entirely in the zero mode sector and is very easily carried out. We must evaluate 

N r oo 

n da m ^ kMmam 

m=l M m =—oo 

(recall M m are the constant units of flux in the U(l) m factors). The summation over M m 
constrain a m to take the discrete values 27 

2nn m 

Oim. — ; • 



26 It is very important below that the ratio of determinants is independent of the fluxes. The following 
comments may point towards the intuitive explanation of this rigorously derived fact. In a background with 
q units of flux, an arbitrary scalar field may be expanded in scalar spherical harmonics which transform in 
the spin q, q+ 1, . . . representations of SU(2). An arbitrary vector field, however, can be expanded in vector 
spherical harmonics which transform in the spin (q — 1, q, . . .) + (q + 1, q + 2 . . .) representations of SU(2). 
The number of vector spherical harmonics with angular momentum q + 1 or greater is exactly twice the 
number of scalar spherical harmonics with the same quantum numbers; as a consequence the contributions 
of these modes cancel in the ratio of determinants. However the number components of in a scalar spherical 
harmonics of angular momentum q is precisely one more than half the number of modes in vector spherical 
harmonics of quantum numbers q — 1 and q respectively. These observations may underlie the reason that 



(4.5) independent of the flux sector. 

We thank O. Aharony and S. Giombi for emphasizing the importance of the summation over flux sectors 



27 

for our problem 
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In other words the integral over «j is replaced by a sum over the uniformly spaced discrete 
values. 

Putting the various pieces together, it follows that the path integral ( |2.5| ) is given by 



4.6 Squashed spheres and genus g surfaces 

Our derivation of the formula (1.6) actually used no feature of the S 2 spatial manifold 



apart from its topology (the topology was used in the assertion that every regular vector 



field can be written in the form (4.4); this is not true when the manifold in question has 
nontrivial one cycles). It follows that ( |1.6j ) applies for an arbitrary metric on the S 2 with 



the volume V2 that appears in (1.7) interpreted as the volume of the manifold. 

Although we will not consider it in the rest of this paper, there is a natural general- 
ization of ( |1.6[) to the partition function of Chern-Simons theory on T> g x S l where £ 5 is 
a genus g manifold of arbitrary metric 28 



N 



rre=l n„ 



n 2sin 



ai(n) - a m (n) \ \ e -V(V) 



(4.6) 



where the summation over n m is restricted so that no two n m are allowed to be equal. A 
straightforward generalization of the derivation of ( |1.6| ) yields (|4.6|) |p2|| . The restriction 
that no two Tim, circ allowed to be equal is presumably because the path integral actually 
vanishes when two eigenvalues of the holonomy matrix are exactly equal. Note that such 
configurations have to be treated specially because our gauge fixing fails precisely where two 
eigenvalues of the holonomy matrix are equal on the sphere. 29 We leave a fuller exposition 
of this point to future work. 

As we have explained in the introduction and will see in detail below, the saddle point 
holonomy distribution of the S 2 partition function is an interesting and nontrivial function 
(with a couple of sharp phase transitions) of v-2 ^ . This interesting behaviour results from 
a competition between potential V{U) that tends to clump eigenvalues, and the measure 
factor in ( jO^ ) which tends to repel them. At higher genus, on the other hand, the measure 
factor in ( |4.6j ) is either absent (for g = 1) or attractive (for g > 1). As a consequence 
the eigenvalue presumably clump at all values of , and the partition function does not 
appear to display any interesting structure as a function of this parameter. In particular 
the high temperature partition function on such manifolds appears always to be given by 
the R 2 partition function regulated by the volume of the manifold under study. 

Unlike the special case of g = 0, the partition function of Chern-Simons matter theories 
on T, g x S 1 for g > 1 has nontrivial dependence on A [21] (and potentially nontrivial 



28 The natural conformal coupling of a minimally coupled scalar to curvature appears to endow these 
scalars with a negative mass for g > 2. This suggests that a theory with a regular scalar is unstable on 
a Riemann surface for g > 2. However the thermal mass appears to stabilize the theory at high enough 
temperatures. 



29 We were not forced to face up to this issue in the special case g = because the measure factor in (f.6) 
automatically kills the contribution of these sectors when g = 0. 
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behaviour) at V2T 2 of order unity. It would be fascinating to find an exact formula that 
generalizes ( |4.6| ) down to temperatures of order unity (of course ( |4.6| ) is only conjectured 
to apply in the high temperature limit, in particular when V2T 2 = O(N)) however we do 
not consider that problem in this paper. 



5. Large iV solutions and level rank duality 

The 'capped matrix model' (1.6) is dominated by a saddle point configuration in the 't Hooft 



large iV limit (N — > 00 with and ^ held fixed). In section |3| above we have determined 
the explicit form of the potential V(U) for several fundamental matter Chern-Simons the- 



ories. In the next section we will turn to a study of the large N limit of (1J3) with V(U) 
given by the specific expressions listed in section ||. In this section, we develop a general 
method to obtain the saddle point eigenvalue distributions for 'integrals' of the form (JO]), 
for an arbitrary single trace potential V(U) (recall the potentials listed in section || are all 
of this single trace form when viewed as a function of the thermal mass and the eigenvalue 
density p; the dependence on thermal mass may be eliminated by extremizing the final 
answer w.r.t. this variable). 30 

5.1 Large N Solution 

As we have explained in the previous section, the summation over allowed values of the 



eigenvalues on in (1.6) is dominated by a saddle point configuration of eigenvalues at large 



N. This saddle point minimizes the following effective potential 

~2 



F([/)-^ln2sin am 2 ai where V(U) = £ V(a m ). (51) 



Of course a m can actually run over only a discrete set of possibilities. As explained in the 
introduction, however, the only effect of this discreteness, in the large N limit, is an upper 
bound on the eigenvalue densities. The eigenvalues within an upper gap saturate a filled 
Fermi sea. Effective Pauli statistics prevents us from infinitesimally varying the positions 
of these eigenvalues. However those eigenvalues that are not located within an upper gap 
have unoccupied slots within which they can move around - effectively continuously in the 
large N limit. The fact that the saddle point extremizes the effective potential (|5.lD then 
imposes the equation 

F'(a m ) = J>t^p ( 5 - 2 ) 
for those a m that do not lie in a filled Fermi sea. In terms of the eigenvalue density function 



p{a) (see (1.8)) we have 

V'(a Q ) = NV j da cot (5.3) 

30 In principle, multi-trace potentials can also be dealt with using the methods of this section, by treating 
the traces that appear on the LHS of the equivalent of ( |5.2| ) below as free parameters, solving the problem 
for arbitrary values of these parameters, and then determining these parameters from the requirement of 
self consistency. We leave a fuller treatment of such potentials to future work. 
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where V represents the principal value (this has its origin in the m / I in (|5.2|)). (|5.2| ) 
applies only for ao in the compliment of the upper and lower gaps of the eigenvalue distri- 
butions. 

In the absence of upper gaps, ( ^.3j ) is a generalization of the Gross- Witten-Wadia 
problem. This problem has been very well studied and admits a general solution in terms 



of an integral involving the potential V(x) (see Appendix B.l for a review). We are 
interested in eigenvalue density functions with upper as well as lower gaps. A simple 
change of variables in ( |5.3D , however, suffices to convert every upper gap into an effective 
lower gap, as we now explain. Let po(a) represent any trial eigenvalue distribution that 
has the right boundary conditions (i.e. is a continuous function that vanishes on all lower 
gaps and equals on all upper gaps). The detailed form of po is arbitrary; we are free 
to choose it to suit our convenience. Under the substitution p{a) = po(a) + tp(a), 
reduces to 



NV J da t/j(a) cot °^ = U(a ) 

da V'(o) = A[p ] 



r /a-e\ (5 ' 4) 

U(a) = V'(a) -NV d6 p (9) cot ( — ^— J 

A[p ] = 1 - J da po(a). 

The important point is that the new variable ip{a) is nonzero only in the compliment of 
both upper gaps as well as the lower gaps (if any). In other words the formulas reviewed 
in Appendix [B.l| - with U and A listed in (|5.4|) - apply to the determination of tp(a). In 



Appendix |B.2| we have explicitly plugged fl5.4| ) into the method of Appendix [B.l| to find 
a formal solution for the eigenvalue distribution of the capped unitary matrix model with 
upper as well as lower gaps. Our final solution is independent of the arbitrary trial density 
po, as of course had to be the case. When V(U) is an arbitrary potential our slightly 
implicit final answer is listed in the appendices. Specifically, the eigenvalue distribution 



is given by ( B.31 ) with i/jy defined in ( B.30 ) in terms of the function &(u)v defined in 



( B.28 ). The contour of integration L ugs runs over all upper gap arcs on the unit circle. 



The locations of these arcs is determined by the normalization condition ( [B.33D . 
5.2 Level Rank duality of the large N solution 

In this subsection we will demonstrate the following. Let p{a) be a solution to the saddle 
point equation 

V'(a ) = NV J da cot "° ~ % (a). (5.5) 



Then 



solves the equation 



p(a) = T ±-(^L- P (a + .)) (5.G) 



V'(a ) = (k — N)V I da cot ^r^p(a) (5.7) 
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where the potential V is defined in ( 2.22Q . In other words the large N saddle point equations 
are consistent with level rank duality. As in the previous section, in this section we assume 
that the potential V(a) takes the single trace form V{a) = Yli^( a i)- However the result 
that our saddle point equations are consistent with level rank duality applies also to multi- 
trace potentials V(U), as we demonstrate in Appendix |D[ 
It follows almost immediately from definitions that 31 



V'(a) 



-V'(a + ir). 



Using ( 5.10| ), ( p.7[ ) (which we wish to prove) may be rewritten as 



-V'(a + ir) 



N 



1-A. 



V I da cot 



ao — a A 



1-A V 2?rA 



1 



p(a + n] 



Using the integral identity 



da cot — ^— — = 0, 



( 5,11 ) reduces to 



-V'(a +tt) = -NV I da cot "° - p(a + Tr). 



(5.10) 



(5.11) 



(5.12) 



(5.13) 



However ( 5.13 ) is simply (f3Jj) evaluated at argument a + it. As we have assumed that 
( |5.5| ) holds at all values of the argument, it follows that ( |5.13 ) and hence (fT?]) also holds. 
On the saddle point, the partition function evaluates to 



- InZ = V[p] - N 2 V / dad/3 p(a)p((3) In 2 sin 



a — f3 



Using the integral identity 



r , , / a-P\ 
/ da In 2 sm 



31 ( 5.10 ) may be demonstrated as follows. Let 

V(U) = A n trU n + c.c. 

n 

so that 



(5.14) 



(5.15) 



(5.8) 



V(a) = J2 A " elUa + c - c 



where A n are arbitrary constants. Now the dual theory potential is given by making the replacement 
tr!7 n -> (-l) n+1 tr!7" in V(U), i.e. 



V(U) = A n (-l) n+1 tr(7 n + c.c 



(5.9) 



so that 



Via) = A n (~l) n+1 e ma =~Y A n e m{a+1T) = -V(a + vr) 



(p.lOj) follows immediately by differentiating both sides of this equation. 
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and the defining relation V[p] = V[p] it is easy to verify that InZ may also be written as 

-InZ = V[p] - (k - N) 2 T J dadPp(a)p{P) In ^2 sin J • (5-16) 

In other words the saddle point equations and the value of the partition function evaluated 
on the saddle point both respect level rank duality. 

6. Exact solution of the S 2 partition function of CS matter theories in the 
low temperature phase 

In the previous section, section ^j, we have recast the problem of computing the S 2 partition 
function of fundamental matter Chern-Simons theories into the evaluation of a 'capped' 
large N matrix model. We have also reduced the problem of evaluating this large N matrix 
model to the determination of solutions of a saddle point equation, and presented a general 
- if slightly implicit- method to determine the solutions of these equations. In order to find 
explicit results for the S 2 free energy of any given matter Chern-Simons theory, we need 
to implement this general procedure for the particular functions v[p] determined in section 
|3[ In this section we partially implement this programme, in a manner we now explain. 

As we have explained in the introduction (and in much greater detail in the context 
of an example in the next section), capped matrix models of the sort that arise from fun- 
damental matter Chern-Simons theories typically have four qualitatively different phases. 
These consist of a low temperature 'no gap' phase, an intermediate temperature 'upper 
gap' or 'lower gap' phase, and a high temperature 'two gap' phase. In this section, we 
present the exact evaluation of the S 2 partition function of all the Chern-Simons matter 
theories studied in section ||] in the low temperature 'no gap' phase. This exercise allows 
us, in particular, to determine the first phase transition temperature of all these theories 
as a function of A. It also allows us to determine the 'critical' value of A, A c for each of 
these theories (the value of A above which the intermediate temperature phase is the upper 
gap rather than the lower gap phase). Finally we are able to verify that the free energy 
of this phase matches smoothly onto the (A independent, non renormalized) spectrum of 
multi traces at low enough temperatures. 

Let the eigenvalue distribution in the no gap phase be given by 



1 1 00 

= 7T + IT E(P» e " <na + P~ne ina )- (6-1) 



2vr 2tt 

n 



The repulsive measure factor in the potential ( |5,1[ ) may be rewritten in terms of p n as 

a i ~ a j at2 f W jo / \ /om ■ a ~ @ 



/7T /*7T 
da / dp p(a)p(P) In 2 sin ■ 
-7T J—TT 



-N 2 



W da f dp P (a)p{P) £ -L + e -***-fl) ( 6 .2) 



XJ 1 



n 
=i 
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where in the last line we have used 

dap(a)e- ina = p n . (6.3) 
It follows that the potential ( |5.1| ) can be written in terms of p as 

OO j 

V s = NT 2 v[p] + N 2 V -PnP-n- (6.4) 

n 

n=l 

As we have seen in section ||, for matter Chern-Simons theories always takes the form 

= - z> — ^ — (P n + ^-«)' 

n=l 

where a n (c, A) are coefficients that depend on a 'thermal mass' constant c which itself is 
determined (by a complicated constitutive equation) as a function of p n . Throughout this 
section we will treat c as a free parameter, solve for p n in terms of c by extremizing the 
potential ( |Q| ) w.r.t. p n to obtain 

Pn = jjra n (c, A). 

(This result is correct only in the low temperature lower gap phase ). We then self con- 
sistently determine c on this solution, using its constitutive equation for c, completing the 
determination of p(a). We now implement this procedure for the various theories we study. 

32 



6.1 Supersymmetric theory 

We first consider the supersymmetric theory. From (3.21) we obtain 
V{p) = V 2 T 2 v[p] 

= — — I c 3 — 6 A / dap(a)Ke [ dyy log tanh Vj^L} a 

6^A V J-n Jc 2 

NV 2 T 2 ( 



dap(a) / dyyl-Y _^_ e -(2m+l)i/ (e i(2m+l)a + e -i(am+l)a 
Jc \ ^0 2m + 1 

? " 3A jH dyy f) ^ T ^ (2m+1)2 '(P(2 m+ i) + P-(2 W+ l))) ) 

(6.5) 



6ttA 

NV 2 T 2 
6vrA 

6vrA l ' 



32 As we have explained in section g, the constitutive equation for c may be obtained by extremizing 
v[p] w.r.t c at arbitrary constant p. Once p n have been obtained by extremization, however, it follows 
immediately that this value of c also extremizes the potential Vs in (3.4). This follows immediately from 
the observation that 

^ = d c V 3 \ Pn +d Pn V 3 \ Pm , c ^, 
ac ac 

together with the fact that d Pn Vs\ Pm ,c vanishes on shell. In other words the value of c may be extremizing 
V, simultaneously V is extremized w.r.t. p n and c. 
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where in the last equation, going from second line to third we have Taylor series expanded 
log terms in y, going from third to fourth line we have used (6.3) and while going from 
fourth to fifth line we have used 

r°° 1 

/ dy ye~ ny = -^e~ nA (l + nA). (6.6) 
J A n 2 



The saddle point for p n ,P-n can be obtained by using ( |BTq ) and extremizing ( |6.4D with 
respect to p- n ,p n . This gives 



P-2n — P2n — 0, 

' Ntt (2n + l) 2 ' 



P-(2n +1 ) = P2n + 1 = ^T^V^^^ 1 + (2n + l)c). 



Plugging ( |6.7| ) into Q6.1| ), we obtain 

m=0 v ; 

Extremizing v[p] w.r.t c at constant p yields ( |3.22D (as explained in the previous subsection, 
this is the same as extremizing ( |6.4| ) with respect to c, on shell) which takes the form 

oo 

5 = A E ^TT e " (2n+1) '^ +1 + 

n= \ oo ( 6 - 9 ) 

W 2 T 2 X 1 2(2n+1)g 

Mr Z^ (2n + i)3 e (l + ^n + ljcj, 

n=0 ' 



where going from first to second line we have used ( p.7|) . 

(|6.9[) completely determine p(a) in the no gap phase. These equations define a legal 
no gap solution provided that everywhere p{a) > and p(a) < s-r- It is easy to convince 
oneself that p(a) is maximum at a = and minimum at a = n. 33 It follows that we have 
a legal no gap solution provided that 

P(0) < -L (6.12) 

and 

P(tt) > 0. (6.13) 



For any c > 0, actually, we can show that 

2 



^ E (^TTF^ 1 + ( 2 ™ + l)5)e- (2 " +1)S < y (6-10) 



because the function is a monotonically decreasing function with respect to c. Hence (6.8) is bounded as 

1 , V2T 2 1 V2T 2 

2^ + W- P(Q) -2^- W (6 ' n) 
So if V 2 T 2 < N, it is legitimate to argue < p(a) < 
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For A < A c it turns out that ( |6.13| ) is violated at a lower temperature than the one which 
breaks ( |6.12| ); the reverse is true for A > A c . One can determine A c for this theory by 
simultaneously solving 



Ktt) = 0, p(0) = — , (6.14) 



which respectively gives 



J_ _ Y^L V i + (2m + i)c (2m+1)g 

2tt vr 2 iV ^ (2m + 1) 2 

m "° (6.15) 
1 , F 2 T 2 ^ l + (2m + l)c „_ f2m+Mg 1 



E 



2vr vr 2 iV ^ (2m + 1) 2 2vrA c ' 

in=0 v 7 

Clearly, A c = \ solves ( |6 . 1 5|) (we could have anticipated this result from the self duality of 
the susy solution). 

In the rest of this section we work on a round sphere of unit radius, so that V 2 = 4-7T. 
For the values A < ^ , the no gap phase transits into a lower gap phase via a Gross Witten 
Wadia type phase transition. We have determined this phase transition temperature 



>Ci 9 W 



v 2 



by numerically solving the equations (|6.9|) and p(ir) = 0. Our result is plotted in Fig |l(a) 
Note that the phase transition temperature increases as a function of A. 

In a similar manner, for values A > ^, the no gap phase transits into an upper gap 
phase via a Gross Witten Wadia type phase at a temperature 



T, = VN 




CugW 



v 2 1 

which we have determined by numerically solving the equations (|6.9| ) together p(0) = 



Our result is plotted in Fig 1(b) . Note that the phase transition temperature decreases 
as a function of A. The phase transition temperature at the critical value of A is T c = 
0.403133 VN. 



It follows from duality that 

N( ug (\) = (k-N)( lg (l-\), 

i.e. that 



C« 9 (A) = ^C/ 5 (1-A). (6.16) 



This relationship is graphically verified in Fig. 1(c) 
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Figure 1: We plot as a function of A for N = 2 supersymmetric case for A < i in Fig. (a), 
A > \ in Fig.(b). In Fig. (a), T denotes the phase transition temperature from no gap phase to lower 
gap phase. In Fig.(b), T denotes the phase transition temperature from no gap phase to upper gap 



phase. We demonstrate (6.16) by Fig.(c 



6.1.1 Low temperature expansion 

On general grounds the free energy of the sphere partition function takes the form 

-lnZ = AT 2 /0^). (6.17) 

A Taylor expansion of the function / may be obtained as follows. To start with, the self 
energy equation (3.22) may easily be solved in a power series expansion in p n . At O(p^), 
c becomes c = (for all A); at higher orders we find 



oo ^ 

2A Yl ^~rj (P(2n+1) + P-(2n+l)) + 0(j%). 
n=0 



(6.18) 



Plugging this solution into ( |6.5[) we find 

NV 2 T 2 °° 



V( P ) 



~ 2^ (2 n + 1)3 ^(2™+!) + P-(2n+l)J 



+ 



2NV 2 T Z X 
~3n 



n=0 
2 \2 



/ oo 1 

L?o(^) (P( 



.19) 



2n+l) + P-(2n+l)) + @(Pn)- 
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The important point in ( 6.1S| ) is that the 0(p} l ) term is independent of A and the 0{p 2 l ) 



term vanishes. Now plugging ( |6.19 ) in (|6.4D and completing square and using the fact that 



(which follow from (6.7)) 



P-2n — P2n — 0, 

V 2 T 2 1 ^/T 6 \ (6-20) 



P-(2n+l) - P(2n+l) ~ ^ (2n + 1)2 + {jjS ) > 



we obtain 



iV y 8 vr V^/ 3 v ^ y V 8 7 V N J \\ N J 

Note that, I ) term is independent of A. Also note that, dependence of A of O 1 



in / (jfi-J is such that under duality, InZ as defined in ( 6.17 ) is invariant. Also note that 
at low temperature, upon setting V2 = 4-7T, O ( j^j term in / (see ( |6.21| )) gives 



^1 

lnZ = ^T 4 C(5), 



where we have used ( |6.17 ). This result is in precise agreement with the high temperature 



limit of the partition function over a gas of (non renormalized) multitrace operators (see 
«). 



6.2 Critical boson 

For the critical boson, using ( 3.12j ) we obtain 



V(p) = V 2 T 2 v[p] = -^-<x 3 + f°° dy / da yp(a) (ln(l - e"^) + ln(l - eT^ )) 



NV 2 T 2 o NV 2 T 2 ^ 



6ir 2ir 
67r 2ir 

"■2 AT\f~Tl r°° 



dy f da yp(a) - ( e - ny+ina + e~ ny - ina ) 
n=i n 



NV0T 2 o NV?T 2 r°° 1 

— s-""— s-y. *»E^fa.+^o 

n=l 

T/9AT 2 , V9AT 2 ^ 1 
= -^7^^ " ^T— E "3 C 1 + + P-n)^, (6-22) 

67T Z7T n J 

n=l 

where going from first line to second, we have used Taylor series expanded log terms in y, 
going from second to third line we have used fl6.3| ) and while going from third to fourth 
line we have used fl6.6|) . 

The saddle point for p n ,p_ n can be obtained by using (|6.22| ) and extremizing (^J) 
with respect to p n ,p- n respectively. This gives 

V9T 2 1 

Pn = P-n = ^~ T ^(l + na). (6.23) 



-40- 



Plugging ( 6.23 ) into Q6.1[) , we obtain 



1 



T 2 V 2 



1 



2vr 2A% 2 ^ n 2 



cos na e 



(6.24) 



n=l 



Extremizing with respect to a at constant /9 n yields the equation ( 3.11| ) we obtain 

oo 1 

a = J2 -e-™ ( Pn + p_ n ) 



n=l 



00 - 



-2ncr 



n=l 



(1 + na) 



(6.25) 



^ 2 T 2 
Nir 



(aU 2 (e-^) + U 3 (e- 2 n) 



where going from first to second line we have used ( gj2f) . The equations (|6T2^ , (|dT24|) ,(|6T2^) 
are independent of A. It follows that the partition function in the no gap phase (and also 
in the lower gap phase) is completely independent of A in this case (see Fig. ( 2(a)| )). 

Note that ( gg§ ) and ( gg| ) are in perfect agreement with Eq.(50) and Eq.(48) of |15| 
provided we choose V 2 = 47r. As in the supersymmetric case, at low values of A the first 
phase transition in this system occurs when p(ir) = 0. This occurs at a temperature 



T c = 0.581067v N. 

As both the no gap phase quantities as well as the condition for the phase transition are 
independent of A, it follows that this phase transition temperature is independent of A. 

Above a critical value of A (denoted by A c ), the first phase transition in this system 
occurs when p(0) = ^x- One can determine A c for this theory by simultaneously solving 



p{0) = 0, P (tt) 



2ir\ r 



and fls!25p . This gives 



A Cri.B. = 0.403033, T c = 0.581067-/W. 



(6.26) 



(6.27) 



For A > A^ ri B - the phase transition temperature depends on A. 

We have numerically solved the relevant equations and plotted the phase transition 
temperature in Fig. |2(a)| and Fig. 2(b). Note that as A is increased beyond A^ ri ' B ', the phase 



transition temperature from no gap to upper gap phase decreases in Fig.2(b) 



6.2.1 Low temperature expansion 

The self energy equation ( p.ll| ) may easily be solved in a power series expansion in p n . At 
C(/9°), a becomes a = (for all A); at higher orders we find 

00 

CT = Y J -(Pn + P-n) + 0(p 2 n ). (6.28) 

71=1 
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Figure 2: We plot phase transition temperature -j= as a function of A for critical bosonic theory 
for A < A° rLB '(= 0.403033) in Fig.(a) and for A > A° ri ' B (= 0.403033) in Fig.(b). In Fig.(a), T 
denotes the phase transition temperature from no gap phase to lower gap phase. In Fig.(b), T 
denotes the phase transition temperature from no gap phase to upper gap phase. 



Plugging this solution into ( |6.22D we find 
NV 2 T 2 



V{ P ) 



°° 1 NV T 2 ( °° 1 \ 

E^(Pn+P-n) + ^~ E" (A. + *>-») +OW- (6-29) 

n=l \n=l / 

As observed in super symmetric case, the 0(pj l ) term is independent of A and the 0(p 2 l ) 
term vanishes. Now plugging ( |6.29 ) in ( J6.-4) ) and completing square and using the fact that 



(which follow from ( 6.23Q ) 



VoT 2 1 /T 6 \ 

we obtain 

^(^)-v(^) 2 « 5 '^(t) 4 (« 3 » 3 (^V + [I^] I- (6.3D 

it / (^fi-J is independen 
regular fermionic theory namely 




Note that / f^H is independent of A. This predicts a simple dependence on A for dual 



Also note that at low temperature, upon setting V2 = 47r, O [jt 1 ^ term in / (see ( 6.31 ) 
and ( |6.17| )) gives 

lnZ = 4T 4 C(5), 

which is exactly same as that for the gas of multitrace operators of free fermions on the 
sphere (recall the critical boson theory has the same single trace spectrum as the free 



fermion theory). For details see ( E.8 ). 

In rest of the section, we simply present final results. All the steps leading to these 
results are completely analogous to those employed in this subsection and the previous one 
(subsections ^T] and |67 
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Regular fcrmion , A <A t 
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Figure 3: We plot phase transition temperature as a function of A for regular fermionic theory 
for A < A^°s- F -(= 0.596967) in Fig.(a) and for A > A^°s F (= 0.596967) in Fig.(b). In Fig.(a), T 
denotes the phase transition temperature from no gap phase to lower gap phase. In Fig.(b), T 
denotes the phase transition temperature from no gap phase to upper gap phase. 



6.3 Regular fermion theory 

For the regular fermionic theory we obtain 

V?NT 2 , V?NT 2 ^ , i 1 

"TT l " A)c 3 - \— ^(-1)™+ 1 -(1 + nc)(p n + p- n )e- nc , 

07TA Z7T z — ' 

V 2 T 2 (-1)™ +1 



V(p) = V 2 T 2 v[p] 



n=l 



P-n = Pn 

p(a) 

(1 - A)c 



, e- nc (l + nc), 
1 VoT 2 v 



2ttN n 

^2 



2vr 2ir 2 N 



£<-« 



1 + m c 



cos ma- 



V 2 T 2 \ yv 1 + m c c _ 2mg 



7riV < — ' m J 



(6.33) 



One can check that under duality last two equations of ( 6.33 ) maps to ( 6.24| ),( |6T25] ) respec- 
tively. The critical value of A and temperature is given by 



A Re g .F. = _ 596967) T = o.477444ViV. 



(6.34) 



The critical value A^ eg ' F ' can be checked to equal 1 — A^ nB ' (see (6.27) for A^? ri ' B ). Upon 



raising the temperature, the no gap solution first transitions to the lower gap solution for 
A < A Bcg ' F ' but to the upper gap solution for A > A Bcg ' F '. Fig. 3(a) demonstrates that the 



phase transition temperature from the no gap phase to the lower gap phase increases as we 
increase A keeping A < A Beg ' F '. On the other hand the phase transition temperature from 

Res' F 

the no gap phase to the upper gap phase decreases upon increasing A when A > A c '. 

For A > A Beg ' F ', the dependence of the phase transition temperature T on A can be 
predicted very simply from duality. Duality maps the regular fermionic theory at A > 
^Reg.F. ^ Q cr j^j ca j bosonic theory at A < A^ n ' B \ In this phase, the phase transition 



temperature, of the critical bosonic theory, is independent of A, as shown in Fig. 2(a) 
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This predicts that phase transition temperature of the regular fermionic theory is given by 



N 



T 1 — X , 

0.581067;/——. (6.35) 



One can check that Fig. 3(b) indeed satisfies this simple relationship. 



6.3.1 Low temperature expansion 



The self energy equation ( |3.7| ) may easily be solved in a power series expansion in p n . At 
O(pn), c becomes c = (for all A); at higher orders we find 

(1 - A)c = A^(-l) n+1 — ( Pn + p. n ) + 0(p 2 n ). (6.36) 

n=l 

Plugging this solution into first equation of ( |6.33 ) we find 



(6.37) 



np) = -^£(-l) n+lJ 3(Pn + P-n) 
n=l 

As the previous case the 0{p] a ) term is independent of A and the 0(p\) term vanishes. 



Now plugging (|6.37) in Q6.4]) and completing square and using the fact that (which follow 



from second equation in ( |6.33D ) 



"-="« = i^(- 1 )" +1 i +0 (^)- < 6 - 38) 



we obtain 



' - S) - - v i (£)'««>4 (3)' <« 3 » 3 <i^ f(£W <-) 



Note that, O term is independent of A and precisely matches with first term in ( 6.31| ) . 

It is easy to verify that, (|6.39|) is consistent with ( |6.32 ). Also note that at low temperature, 
setting V2 = 47r, we obtain InZ = 4T 4 £(5), which is exactly same as that for the gas of 



multitrace operators of free fermions on the sphere. For details see ( |E.8| 
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6.4 Regular boson 

For regular bosonic theory we obtain 




VT 2 1 

^2 



1 T V \— - 1 

h — tt-4 ^ cos(na)e" n<J (l + no) 



2vr 2Ntt 2 

n=l 

n=l 

(6.40) 

For this model, the critical value of A, which is denoted by A c ', depends on A6- When 
A6 — > oo, results are identical to critical bosonic theory. A few other values of A c ' are 

A Re g .B. = 0.360884, T c = 0.555125\/iV, A 6 = tt 2 . 

A Re g .B. = .373102, t c = 0.559636ViV, A 6 = 8^ 2 . (6.41) 
A Re g .B. = 0.382887, T c = 0.564807\/iV, A 6 = 32vr 2 . 

Res B 

As in previous subsections, for A < A c ' the no gap solution first transits to the lower 



gap solution (see Fig. 4(a) ); the phase transition temperature increases as a function of A. 



R,Gf£ B 

On the other hand for A > A c ' the no gap solution first transits into the upper gap 
solution (upon raising the temperature); the phase transition temperature decreases as A 
is increased (see Fig. 4(b)). 



6.4.1 Low temperature expansion 



The self energy equation (|3.4| ) may easily be solved in a power series expansion in p n . At 
O(pn), o becomes o = (for all A); at higher orders we find 



1 + 



\ oo 

\o = J2-(Pn + P-n)+0(p 2 n ). (6.42) 



Plugging this solution into first equation of (|6.40| ) we find 



v( > NV 2 T 2 ^ 1 



2tt ^— ' w 

n=l 



+ iVV^ 



12vr fY 6 



5.43) 
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Regular Boson at Ag=7T ,A>A C 

Regular Boson at A 6 =n 2 , A <A ( . 7 




0.5530 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 A - 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 <- A 

0.0 0.1 0.2 0.3 0.4 0.4 0.5 0.6 0.7 0.8 

(a) (b) 

Figure 4: We plot the phase transition temperature -7= as a function of A for regular bosonic 
theory for A < A^ eg B (= 0.360884) in Fig.(a) and for A > A^ egB (= 0.360884) in Fig.(b). In Fig.(a), 
T denotes the phase transition temperature from no gap phase to lower gap phase. In Fig.(b), T 
denotes the phase transition temperature from no gap phase to upper gap phase. 



As observed in the previous models the 0{p] l ) term is independent of A and the 0{p 2 l ) 
term vanishes. Now plugging ( |6.43 ) into (|6.4j ) and completing square and using the fact 
that (which follow from second equation in Q6,40D ) 

- - ^WofSV (,44) 



-n — Fn 



2NtTT1 2 \N 3 

we obtain 

> m - v i & -4 m -4 ♦ ^ 2 5 

(6.45) 

Note that the O (^j^j term is independent of A and infact agrees perfectly with the gas of 
multitraces of the free boson, flE.4|) . 

6.5 Critical fermion 

For the critical fermionic theory we obtain 

o V2NT 2 o VoNT 2 A ±1 1 

V(p) = V 2 T 2 v[p] = (1 - A + g)c 3 ~ E(" X ) ^ C 1 + n£ )(Pn + P-n)e~ nc , 

n=l 

p. n = Pn = 2 { > e~" c (l + nc), 
zirJy n z 

1 V2T 2 ^— k. 1 + m c _ - 
p(a) = 7T— > (-1) cos ma = — e , 

m=l 



V2T 2 \ V 1 • n, r .,„ 
ttN 



m=l 

g(\) = j. (6.46) 
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Critical fermion , A >A t . 

Critical fermion ,A<A ( . T 




0.38 1 X — ' 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 A 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.6 0.7 0.8 0.9 1.0 

(a) (b) 

Figure 5: We plot the phase transition temperature as a function of A for critical fermionic for 
A < A° rlF (= 0.639116) in Fig.(a) and for A > A^ rLF '(= 0.639116) in Fig.(b). In Fig.(a), T denotes 
the phase transition temperature from no gap phase to lower gap phase. In Fig.(b), T denotes the 
phase transition temperature from no gap phase to upper gap phase. 



One can check that under duality third and fourth equation of ( |6.46| ) maps to third and 
forth equation of ( |6.40| ) with the identification g(\f) = , 2 ^ b One can check that, given 



the A&, Ag in ( |6.41| ), the A)? riF ' for the critical fermionic theory satisfies A^ n,F ' = 1 — A^" eg B '. 
For example for A* eg,R = 0.3731, Ag" = 8vr 2 one gets g = 0.69913 and we obtain A^ ri F - = 
0.6269. For A < \° liF - the no gap solution first transits to the lower gap solution upon 
increasing the temperature; the phase transition temperature increases as A is increased 
(see Fig. 5(a)| ). On the other hand for A < A^ ri F - the no gap solution first transits to 



the upper gap solution upon raising the temperature; the phase transition temperature 
decreases as A is increased (see Fig. 5(b)). 



6.5.1 Low temperature expansion 

The self energy equation ( f3.18|) may easily be solved in a power series expansion in p n . At 
O(pn), c becomes c = (for all A); at higher orders we find 

oo 1 

(l-\ + g)c = \ ]T(-ir +1 - ( Pn + p- n ) + 0{pl). (6.47) 

n=l 

Plugging this solution into first equation of fl6.46 ) we find 



n=l 

nv 2 t 2 x 2 AA, ^^i . A 3 _ 4 

I2J- 1 ) n+1 -{Pn + P-n)\ +0{p 4 r , 



(6.48) 



+ 



12tt(1-A + ^ 



As before the 0(pl l ) term is independent of A and the 0(p 2 1 ) term vanishes. Now plugging 
(6.48) in (|6.4[ ) and completing square and using the fact that (which follow from second 
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equation in (|6.46|) ) 



V 2 T 2 




(6.49) 



P-n = Pn = 



2Ntt 



we obtain 




perature limit of the gas of multitrace operators of the free boson theory (E.4). 



One can check that the O ((^r) 4 ) term in ( |6.50D is related to the term of the same 



order in ( |6.45 ) by duality. 



7. Exact solution of the Large N capped GWW model 

In the previous section we have solved for the eigenvalue distribution p(a), and thereby 
determined the free energy of systems we have studied, in the low temperature no-gap 
phase. It is much more difficult to solve the saddle point equations in the intermediate 
temperature lower or upper gap phases, and the high temperature two gap phase. In this 
section we accomplish this task; not for the realistic potentials v[p] described above, but 
instead for a toy model potential 



In this case V(a) = —N( cosa. 

Our toy model shares the following key qualitative feature with the function v[p] that 
comes from real matter CS theories; the potential V[U] always has a minimum at U = I, 
and the depth of the potential increases as a function of £. We hope that our solution of 
one and two gap phases for our toy model, presented in this section, will aid in the solution 
of the saddle point equations in the same phases for the potentials v[p] that arise in real 
matter CS theories; however we leave the investigation of this issue to the future. 

7.1 No gap solution 

The 'no gap' solution of these equations is identical to the same solution in the uncapped 
GWW model and famously takes the form 

1 + Ccosa 

= ^ • ( ' 

For all A < 1, the no gap solution is the correct saddle point for small enough £. The no gap 
solution ceases to be a saddle point of our model if either p(a) goes negative somewhere or 
if p(a) exceeds somewhere. Which of these occurs first (upon monotonically increasing 
A) depends on the value of A. 
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p(a) takes its maximum/minimum value at a = O/ir; these values are given by 
It follows that if A < \ the no gap solution goes negative before it hits the upper bound. 
In that case the intermediate ( phase has a single lower gap. On the other hand if A > ^ 
then the intermediate C, phase has a single upper gap. 

7.2 Single lower gap solution 

This is the intermediate ( phase that occurs if A < 5. This phase kicks in for ( > 1. The 
eigenvalue distribution in this phase is same as that for the standard (uncapped) GWW 



solution (see Appendix B.1.2 for a review) and is given by 



Ccos(§) fl . 2 a . 2 a 1 

p( a ) = sin - for sin - < -, 

7T V C 2 2 Q 

p(a) = for sin 2 — > — . 
HK ' 2 C 

Note that we always have a lower gap (a region where the eigenvalue distribution vanishes) 
as C > 1- The eigenvalue distribution is maximum at a = 0, and the maximum value is 
given by It exceed the maximum permissible value at C = > 1 (recall 2A < 1). At 
this point the system undergoes a second phase transition to the one lower gap and one 
upper gap phase. 

7.3 Single upper gap solution 

This is the intermediate C phase when A > ^. This phase kicks in at C = t — 1 < 1- I n 
Appendix |C.1| we have solved for the eigenvalue distribution in this phase using the general 



analysis of Appendix R^. In this subsubsection we give a much simpler derivation of the 



eigenvalue distribution in this simple phase; of course our final result agrees with that of 



Appendix C.l 



In the case of phases with only one upper gap the general analysis of capped matrix 
models presented in Appendix |B.2| greatly simplifies if we make the following natural choice 
for the trial eigenvalue distribution function po(ce)', 

1 



2ir\ 

As the integral of the cotangent against a constant vanishes, with this choice we find that 
U(a) = V'(a) in ( ^.4[ ). It follows that the integral equation in ( |5~i| ) reduces precisely 
to the Gross Witten integral equation. The boundary conditions on tp are the converse 
of the GWW boundary conditions; ip{a) = for |a| < a. We can turn these boundary 
conditions into the GWW boundary conditions by working with a new variable 9 = a — tt. 
We have ip(0) = for |#| > ir — a. When rewritten in terms of 9 the GWW force function 
U{a) = NC,sma turns into — NQsmO. 

It follows that —ift(9) is exactly the GWW eigenvalue distribution for the potential 
V(a) = N^cosO with one minor difference; the eigenvalue distribution density is normal- 
ized so that „ 

(-ip)da = --1, 
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(the RHS is unity in the standard generalized GWW problem). This change is easily 
accounted for (see Appendix || below) and yields 



7T 



Substituting 9 = a — n above we have 



1 Jsin(f)| .a . 2 a £-1 

= ^ fOT C ° S 2 > ~C~- 



(7.3) 



This solution has an upper gap region as we have assumed that £ > v — 1. The minimum 
value of the eigenvalue density occurs at a = it and is given by 



1 / 1 

2tt 



The minimum goes to zero for 

at which point the system transits to the one upper gap and one lower gap phase. 
7.4 One lower gap and one upper gap solution 



In Appendix |C.2| we have demonstrated that our toy model also admits a solution with one 
lower gap and one upper gap. The upper gap extends on the arc on the unit circle from 
e~ ta to e m and the lower gap extends on the arc on the unit circle from e lb to e~ lb (see 
Fig. 9(c) in the Appendix) where a and b are positive and are determined by the equations 



4vrA / „ / . 2 a -2 a 1-2 b -2 a ^' 



sin 2 — sin j y sin | — sin 2 



1 l' a r cos a C , 

da — — j^^^=^^^= = 1 H — (cos a + cos b) 



(7.5) 



47rA / „ /■ 2a ■ 2 a / ■ 2b • 2 a 2 

• y_a wsui | - sin f ysin | - sm | 
In the compliment of these two gaps the eigenvalue distribution is given by 



| sinal // . 2 a . 2 a\ ( . 2 6 . 2 a\ , 

47T/9(a) = — 4 / I sm sm — I I sin sm — l\ where 

7r A V \ 2 2/V 2 2 



I 



a 



dQ 



(7.6) 



!ciis« - cos n ) x (sin 2 | - sin 2 |) (sin 2 | - sin 2 £ 



By expanding ( C.16| ) in Appendix C.2 near a = b we find (see (|O30|) in Appendix |03l|) 



p(a) = Byfb--a + 0(6 — a) 2 where -B is a known positive constant. Similarly, near a = a 
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p(a) = + D\Ja — a + 0(a — a) 2 where D is another known negative constant (see 
QC.37| ) in Appendix C.3.4). In other words the eigenvalue density in the 'cuts' continuously 



matches onto the eigenvalue densities in the upper and lower gap regions. 

Over what range of £ does the double cut solution exist? In the double cut solution, 
a and b are functions of £ and A. For some formal purposes it proves useful to invert this 
relationship and regard £ and A as functions of a and b. As analyzed in the Appendix |C.5| , 
C is a decreasing function of b at fixed a while A is an increasing function of b at fixed a. 
By numerical evaluations, we can also guess that both £ and A are increasing functions of 
a at fixed b. From these facts, we can obtain lower bounds on the range of £ for which the 
double cut solution exists by setting a = and b = it. 

Let us first set b = it in which case (17. ST) reduce to 



a 



£ = / da 



4irX 7 -a ^(sin 2 f - sin 2 f )(1 - sin 2 f ) 2Xc ™l 

1 c , 1X 1 r , cos « 1 ( , 1 

1 H — (cos a — 1) = — - / da — = — 11 

47rA - / - /(sin 2 f -sin 2 f)(l-sin 2 f) A V 2co Sl 



(7.7) 



( |7.7j) are easily be solved to get determine a and £ in terms of A; we find 

C= 1 

4A(1-A)' (7g) 
cos^ = 2(1 - A). 

Note that ( |7.8| ) matches exactly with the limit of single upper gap solution at the edge of 
its validity (see (ffl|)). 

Let us next consider the limit a — > 0. In this limit the integration range in ( [T.5| ) goes 
to zero, but the integrand simultaneously diverges in such a way that the integral is finite. 
The equations reduce to 

c 1 



This implies that 



2Asin|' 

2 (7.9) 

l + |(l + cos6) = — — 
^ 2A sm I 



C = ^, siu^ = 2A, (7.10) 



which again matches with the limit of single lower gap solution at the edge of its validity. 

It follows that the double gap solution studied in this section exists if and only if 
C > ^ and C > 

4A(i-a) ( a l so see Appendix |C.5|) . When A < \ the first inequality is more 
stringent, and the double gap solution exists only when ( > i.e. only when the GWW 
single lower gap solution ceases to remain a valid solution. On the other hand when A > \ 
the second inequality is more stringent, and the double cut solution exists only for £ > j™, 
i.e. only when the single upper gap solution ceases to be valid. 
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The discussion above makes it plausible that when A > | the double gap solution of 
this subsection becomes identical to the single upper gap solution when £ = • ^ 

also suggests that when A < g, the solution of this subsection becomes identical with the 
single lower gap solution at ( = j^. In Appendix C.3.2 and Appendix C.3.1 we have 



demonstrated that this is indeed the case by explicitly computing the double gap eigen- 
value distribution in these limits, and comparing with the respective single gap eigenvalue 
distributions. 



To end this subsection we consider the limit of large £ at fixed A. In Appendix |C.3.3 
we have demonstrated that in this limit 

a = n\ - 4sin(^A)e-^ Asin(7rA) + <D[( ^CAsin^A)^ 

b = vrA + 4sin(7rA)e-^ Asin(7rA) + 0[( e -^M^\ 

1 _ x /a -a 
p(a) = —5— cos 



7r 2 A V b — a 

Note in particular that the eigenvalue distribution tends exponentially to the distribution 



( 1.10 ) in the limit of large £. 



7.5 Level Rank Duality 

Recall that 

V(TrU n ) = V ((-l) n+1 TvU n ) . 
For the GWW toy model of this section 

V(U) = V{U) = -N^ (TrU + TrU^ =-(k- N) 2 ^ ^ (TrU + TrU^ . (7.12) 

In other words the GWW model is self dual under level rank duality together with the 
appropriate rescaling of £. Earlier in this section we have explicitly determined the saddle 
point eigenvalue distribution for the matrix model with the potential V(U). This eigenvalue 
distribution is a function of A, £ and of course a (the phase on the unit circle). Let us 
denote our saddle point eigenvalue distribution by p(\,(,a). Level rank duality predicts 
that the exact solution to this model, obtained in the previous section, obeys the identities 

,(<! - X),^,a) - ^ (-L -p(A, C , Q + ,)) . (7.13) 



In Appendix |C.4| we have verified in great detail that ( 7.13j ) is indeed correct 



7.6 Summary 

We end this section with a summary of the behaviour of the capped Gross- Witten-Wadia 



model. Fig. 6(a) summarizes different phases of the GWW model. 
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7.6.1 A < \ 



For £ < 1, the eigenvalue distribution is wavy (i.e. has no gaps) and is given by 
Fig-gi 



for a plot at A = 4 and £ 



For 1 > C > 



□]) (see 

r , the eigenvalue distribution has 
a lower gap centered about tt and is given by ( |7.2[) (see Fig. 7(b)| for a plot at A = 4 and 
£ = 2). For £ > ttt, the eigenvalue distribution is given by (|7.6|) which has a lower gap 



l 

4A 2 



centered about it and an upper gap centered around zero (see Fig. 7(c) for a plot at A = \ 
and £ = 4.6). At large values of (, the eigenvalue distribution tends to ( l.lCf ) (see ( 7.11| ) 
for more details, and see Fig. 7(d) for a plot at A = \ and Q = 11.03). 



7.6.2 A > \ 

i 



For £ < j — 1, the eigenvalue distribution is wavy (i.e. has no gaps) and is given by ([7. If ) 



(see Fig. 8(a) for a plot at A = 0.51 and £ 
>er gar 

0.51 and ( = 0. 



Fori 



1 > C> 



i 



-, the eigenvalue 



4A(1-A) ' 

distribution has an upper gap (or plateau) centered around and is given by ( |7.3| ) (see Fig. 

1 the eigenvalue distribution is 



8(bl for a plot at A = 0.51 and ( = 0.98). For C > 4A(l _ A) , 
given by ( |7.6| ) which has a lower gap centered about 7r and an upper gap centered around 
zero (see Fig. |3(c) for a plot at A = 0.51 and ( = 1.13). At large values of £, the eigenvalue 
distribution tends to ( 1.10D (see ( |7.11| ) for more details, and see Fig. S(d) for a plot at 
A = 0.51 and ( = 5.06). 



8. Discussion 



In this paper we have presented a wealth of exact results for the thermal partition function 
of several matter Chern-Simons theories as a function of temperature and A in the strict 
large N limit. Our results obey some non trivial consistency checks; they are in perfect 
agreement with several conjectured Giveon-Kutasov type dualities between pairs of theo- 
ries. In every case the partition function also reduces, at low temperatures, to the N and 
A independent partition function of a gas of multi trace operators, in agreement with the 
nonrenormalization theorem for multi trace operators in these theories) 1£, 18]. 
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(a) 



(b) 



(c) 



(d) 



Figure 7: The eigenvalue distribution of the capped GWW model at A = 0.25 for increasing values 
of (: Fig.fT(a)1 is for C = \, Fig ff(b)| is for ( = 2, Fig.|7(c)|is for C = 4.6 and Fig.f^djs for C = f f .03. 
In this case as we increase C, eigenvalue density distribution first develops a lower gap as Fig. [7(b)| . 
Further increasing eigenvalue density distribution develops a upper gap as well Fig. 7(c). 



While the consistency checks reported in the previous paragraphs give us confidence in 
our results, some aspects of the procedure employed in this paper require fuller discussion. 
To start with, it would be useful to have a more thorough justification for using ( |1.2|) as the 
starting point of our analysis. Subleading corrections (in ^) to (F2) could modify ( |1.(]| ) 
in several ways; 34 for instance they could spread the discrete delta functions in ([O]) into 
a distribution function for eigenvalues that is peaked about «j = with a width 5 that 
goes to zero in the large N limit. With such a spread out eigenvalue distribution function 
two eigenvalues can presumably occupy the same slot. 35 Hamiltonian methods indicate 
that the contribution of configurations with two nearly equal eigenvalues is very small. 
Consequently the path integral must include a mechanism that suppresses contributions 
of configurations with nearly equal eigenvalues; it is important to understand in detail 
how this works. Relatedly, it would be useful to thoroughly understand the reason for the 
exclusion of configurations with equal eigenvalues in the path integral (^1]) on S g x S 1 
where £ s is a genus g surface. We hope to return to these important issues in the future. 

The procedure we have employed in the computations of this paper is a mix of two 
different methods. We used lightcone gauge on R 2 to compute the effective potential v[p], 
and then proceeded to evaluate the Chern-Simons path integral supplemented with this 
effective action in temporal gauge. It would of course be more satisfactory to rederive 



34 



We thank O. Aharony and J. Maldacena for discussions on this point. 
35 The energy price from the Vandermonde in the measure, -In S, is presumably negligible compared to 
the energy gain from two eigenvalues occupying the same slot. 
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(a) 



(b) 



(c) 



(d) 



Figure 8: The eigenvalue distribution of the capped GWW model at A = 0.51 for increasing values 
of (: Figfgjs for C = |, Fig.gbjs for C = 0.98, Fig.gcjs for C = 1.13 and Fig.|(djs for ( = 5.06. 
In this case as we increase £, eigenvalue density distribution first develops a upper gap as Fig.|S(b)|. 



Further increasing eigenvalue density distribution develops a lower gap as well Fig. 8(c) 



all our conclusions from computations that utilize temporal gauge all the way through. 
This would also allow the evaluation of the partition function of Chern-Simons theory on 
arbitrary genus g surfaces at finite values of the temperature] 21, 25 1. The main obstruction 
to this procedure are a set of spurious divergences that appear in computations in matter 
Chern-Simons theories in temporal gauge [16]; it may prove possible to understand these 
divergences well enough to deal with them. We hope to return to this point in the future. 

Putting aside all these concerns, the work presented in our paper raises several inter- 
esting questions. One of the most straightforward of these relates to the possible existence 
of a wider class of dualities than those verified in this paper. It should be straight forward 
to generalize the computations of v (U) presented in this paper to a wider class of matter 
Chern-Simons theories and perhaps discover new dualities. 

Giveon Kutasov dualities continue to hold away from the strict vectors like limit; 
for instance they continue to hold at large N c and Nf with held fixed. It would be 
interesting to investigate whether the same is true of nonsupersymmetric dualities, perhaps 
by performing a computation of the S 2 partition function at first order in j^-. 

Relatedly, in this paper we have focused on theories with a single Chern-Simons gauge 
group and matter in the fundamental representation. The story is even richer in Chern- 
Simons theories, like ABJ theory, with two gauge groups - U(N) and U(M) and bifunda- 
mental matter. At least in the limit in which ^ is small the eigenvalues of U(N) should 
develop upper and lower gaps (see [17] for an analysis of the partition function of the free 
theory). It would be very interesting to analyse how the situation evolves as a function 
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of as the corresponding Vasiliev theory turns into a vanilla two derivative theory of 
gravity [L7| (IIA theory on AdS^ x CP 3 with 2-form flux) at large A when ^ = 1. 



It would be interesting, and should be possible to generalize the computations presented 
in this paper to the computation of the partition function of the relevant theories on S 3 . 
The results of such a computation must agree with the results obtained through localization 



for supersymmetric theories [42|, but will, of course, be more general than those results. 

In section ^, we have computed eigenvalue density distribution in the no gap phase for 
several Chern-Simons matter theories. Following section [7|, it would be very interesting to 
find out eigenvalue density distribution in the single lower gap, single upper gap and one 
lower-one upper gap phases of these theories. 

An interesting aspect of the S 2 partition functions determined in this paper is the 
presence of phases with upper gaps. Such phases are absent in the partition function of 
Yang Mills theories As we have seen below, the presence of upper gap phases is a 

direct consequence of the discretization of the eigenvalues of the holonomy matrix, which, 
in turn, is a direct consequence of the summation over flux sectors in the path integral. 
What is interpretation of this new phase in terms of states on S 2 ? As the existence of this 
new phase relies on the summation over fluxes, it is natural to guess that the new phase is 
one in which states with fluxes on the S 2 (states dual to monopole operators) proliferate in 
the system. Note that monopole operators in Chern-Simons theories all have dimensions of 
order k which is very large in the large TV" limit. Such operators have essentially no impact 
on the partition function at temperatures of order unity. However such operators could 
proliferate at VT 2 = £iV as studied in this paper; our results indicate this happens when 
an upper gap is formed, i.e. at £ ~ at small A. It would be interesting check whether 
or not this is true, and to flesh it out. 36 

It has recently been conjectured that level k Chern-Simons matter theories with mat- 
ter in the fundamental / bifundamental representations, admit a dual description governed 



by parity violating Vasiliev's higher spin equations [43, |44[] (See |45| for more detail and 
references) at finite values of the 't Hooft coupling A = ^ |l6|, [T^j . The thermal system 
studied in this paper is thus dual to the Euclidean Vasiliev system in global AdS space 
compactified on a circle on circumference 2ttR = The field theory analysis performed 
in this paper implies that this Vasiliev system admits a classical limit in large N limit 
with V2T 2 = £iV and ( held fixed. The saddle points obtained in this paper should map to 
classical solutions of this bulk description; however the equations of motion governing this 
classical description are not necessarily Vasiliev's equations. This is because quantum cor- 
rections to Vasiliev's equations proportional are not suppressed compared to classical 
effects in the combined high temperature and large N limit studied in this paper; such 
terms could modify Vasiliev's classical equations. It would be fascinating to determine the 
effective 3 dimensional bulk equations dual to the large N limit of this paper, and to study 
the classical solutions dual to the saddle points described in this paper. 

Finally, one of the most interesting aspects of the study of Chern-Simons matter the- 
ories is that we have tantalizing hints that miraculous but well known properties of super- 



°We thank S. Shenker for discussions on this topic. 
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symmetric theories - level rank type dualities - are also true of nonsupersymmetric theories. 
There is a more interesting, more miraculous and better known property of the supersym- 
metric ABJ theory, namely that it reduces to a theory of gravity at M = N and A — > oo. A 
really interesting question is whether similar result holds for a non supersymmetric theory, 
for instance the theory of minimally coupled bifundamental fermions. The study of this 
question sounds like a fascinating programme for the future. 
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A. Interpretation of the Fadeev-Popov Determinant 

In this Appendix we present a heuristic demonstration that this FP determinant simply 
converts the path integral over the two dimensional hermitian zero mode field a(x) into a 
path integral over the two dimensional unitary field U(x). The operator in the determinant 
has no x derivative, so the determinant is a product of determinants, one for every x (this 
is the part of the evaluation that is heuristic - strictly we need to regulate to make sense 
of this statement). Specializing for a moment to the group U(N), the determinant is given 



The product is taken over all (positive and negative) integer values of n when m / I but 
only over all nonzero integer values of n when m = I. Actually all terms in the product 
are independent of a (and so can be absorbed into the normalization of the path integral) 
when m = I. The nontrivial contribution is from m 7^ I. 
After removing an a independent factor we have 



where we dropped overall constant regulated suitably. The last product vanishes at a m — 
ai = 2irn for every integer n. It is thus plausible (and true) that it is proportional to 



by 
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Consequently, the Fadeev-Popov Determinant is heuristically given by 

a m - ai 



Det' s (d 3 + [A 3 ,}) = nil 2 sin 

x m^l 



nn 

x m<l 



where we put the factor 2 so that (see §5.2) 

J DA 3 Det' s (d 3 + [As,]) = J DU(x), 

where U(x) = e^ As and DU is the path integral over U(x) with the left-right invariant 
Haar measure over U(x). The Haar measure is consistent with level rank duality. 



B. Solution to the saddle point equations of 'capped' unitary matrix mod- 
els 

As we have explained in the main text, Chern-Simons theories coupled to fundamental 
matter fields in the 't Hooft limit are governed by the saddle point equations of a large 
N unitary matrix integral with one additional restriction; the density of eigenvalues is 
bounded from above at ^- r. In this appendix we present a general method to solve these 



'capped' large N unitary matrix integrals. In subsection B.l we first review the general 



method to solve ordinary (uncapped) unitary matrix models following the method of [11]. 



In subsection |B.2| we then proceed to use a very similar method obtain an exact solution 
to the saddle point of capped Unitary matrix models. The formulas presented in this 
appendix apply, with great generality, to any single trace capped unitary matrix integral; 
however the final results are a little implicit. In appendix ^ below we show how this formal 
solution can be converted into explicit formulas in the context of a simple example. 

B.l Review of standard Unitary matrix integrals at large N 

As we have explained in the main text, the large TV" eigenvalue density function of the 
standard (i.e. uncapped) unitary matrix integral 



/ DUe- v M 

obeys the integral equation 



/ 



NT / da ^(a)cot ( — — - ) = U(a ) 

J \ 2 J ( B1 ) 

da ip(a) = A 

with U(a) = V'(a). In this subsection we review the standard solution of the integral 



equation ( B.l ). In the next subsection we will use this solution to obtain the solution to 



the capped unitary matrix integral. 
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B.l.l General solution of (gl|) 



In order to obtain the general solution of (|B.1| ) we find it useful to use the following complex 
variables (see for eg. fll|1 ) 



Note that 



z = e ia °, A i = e ia ', Bi = e 



dz ao — a zq + z 
da = — , cot = i- 



ib, 



iz 2 zq — z 

(B.l) may be rewritten in these variables as 



NV I °-—^(z) = U(z ), 

- -u - (]}.2) 



/ %{z) = A 
J iz 



where the integrals in flB.2|) run clockwise over the unit circle in the complex plane. 

Let us suppose that the solution to this equation, ip(z), has support on n connected 
arcs on the unit circle in the complex plane. We denote the beginning and endpoints of 
these arcs by Ai = e lCLi and Bi = e lhi (i = 1 . . . n). Our convention is that the points A\, 
B\, A2, B2 ... A n , B n sequentially follow each other counterclockwise on the unit circle. 
We refer to the n arcs {Ai, Bi) as 'cuts'. The arcs (Bi, Ai + \) (as also (B n , A\) ) are referred 
to as gaps. By definition, tp has support only along the cuts on the unit circle. 

We use the (as yet unknown) function ^(z) to define an analytic function on the 
complex plane 

'~ B * dzu + z 
'A, 



, / x f Bi dzu + z , . . _ . 



where the integral is taken counterclockwise along the n cuts, (Ai, Bi) of the unit circle in 
the complex plane. It is easy to see that the analytic function <3?(u) is discontinuous along 
these n cuts. Let $(z) + denote the limit of $(«) as it approaches a cut from \u\ > 1, and 
let &~(z) denote the limit of <I>(n) as it approaches a cut from |n| < 1. 37 Then 

<S> + (z) - $-(z) = 4inp(z). (B.4) 



37 We use analogous notation for other analytic functions below. Completely explicitly, let z denote a 
complex number of unit norm. The symbol F + (z) will denote the limit of F(u) as u — > z from above (i.e. 
from |m| > 1), while F~(z) is the limit of the same function as u — > z from below (i.e. from \u\ < 1). Note 
that along a cut F~(z) = —F + (z). 
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38 On the other hand 

" Bl duj z + u , , W(z 

■ / ■ I ~ = — HTF 



<&+(*) + *"(*) = 2V> / ^1^^) = ^^ (B.6) 



(the first equality is the definition of the principal value, while the second equality follows 



using (|B.2| )). Moreover it follows immediately from (£^2) and (|B.3| ) that 

lim $(«) = A. (B.7) 



As the RHS of ( |B.6j ) is a known function, we are posed with the problem of determining 
<3?(u) given its principal value along a cut. The problem of determining an analytic function 
from its discontinuity along a cut is, of course, standard. A simple variable change allows 
us to turn the principal value of $(it) into the discontinuity of a new function, as we now 
explain. 

Consider the function 



h(u) = s/{Ai - u)(B 1 - u){A 2 - u)(B 2 - u) . . . {A n - u){B n - u). (B.8) 

We define h(u) to have cuts precisely on the n arcs on the unit circle that extend from 
(Ai, Bi). This definition fixes the function h(u) up to an overall sign. This sign will cancel 
out in our solution for <£> below, and so is uninteresting. For future use we note that when 

u = e ia 

h (n) = [[ 4e 2 e I sm I J - sin I — J J . (B.9) 

39 We use the function h(z) to define a new function, H(z), via the equation 

= h(z)H(z). 

Using the fact that h + (z) = —h~(z) along the cut, QB,6|) turns into 



,i8 (B.4) may be derived as follows. Let D(u) represent any function that is analytic in a region that 
encloses all n cuts. The integral 

' du<f>(u)D(u) (B.5) 



over a contour that encloses all n cuts is, by definition given by J dzD(z)(& + (z) — <fr~ (z)) where the 
integration runs anticlockwise along the n cuts. On the other hand, by substituting (B.3) into (B.5), 
interchanging the order of integration, and evaluating the contour integral over u by use of Cauchy's 
theorem, we find the alternate expression J dz4:irip(z)D(z). The equality of these two expressions for 
arbitrary D(z) implies (B.4). 
39 This may be seen by using 

/ iai ioy % iav „ ( ( «1 + i>l \ ( <l\ — b\\\ 

(e 1 — e ) (e 1 — e ) = 2e V 2 / I cos I a I — cos I I I 

4e<^ + «) fsin 2 ( a -±-±) - sin 2 f f - ' ^ 
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Note that at large v, H(v) = O(-n) so that 



Co 



dv H ^ . = (B.12) 
2m(v -u) y 1 



where the the contour runs counterclockwise over a very large circle at infinity. As- 
suming that the function H(u) has singularities only along the cuts (Ai,Bi) and that the 
function U(u) is analytic on the cuts, the use of Cauchy's theorem on ( B.12Q yields 



Jc cut3 2ni(v-u) nJ Larca Nh+(z)(z-u) 2tt J CcuU Nh(v)(v - u) 

(B.13) 

In the equation above, v is a variable on the complex plane while z is a variable on the unit 
circle of the complex plane. The contour C cu t s encloses each of the n cuts (Aj, -Bj) but does 
not enclose the point u. The integration region Lares runs counterclockwise along n cuts 
on the unit circle, integration region for each cut is from Ai to Bi. The second equality in 
( B.13| ) is obtained by using ( |B.11| ). The third equality uses h + (z) = —h~(z) together with 



the assumption that U(z) has no singularities on the n cuts. 

In practical situations, ( |B.13|) may itself be evaluated using contour techniques. Sup- 
pose, for instance, that U(z) is a meromorphic function whose poles, located at the points 
z\ . . .z r . Then we find 

jj / \ If, U(z) iU{u) ^ U(z) 

H(u) = — / dz — — — t-J— -V4 - > iReSz=z h — r-r-rr 1 r ( B -14) 

y> 2vr J c , Nh(z)(z-u) Nh(u) ^ k Nh(z)(z-u) y ' 

where the contour C runs counterclockwise over an infinitely large circle. In simple exam- 
ples the last contour integral either vanishes or evaluates to something simple. 

B.1.2 The GWW problem 

As an example let us consider the GWW matrix model[10, 12, 11]. Let the potential be 
given by 

V(U) = -N^ (TrU + TrC/ 1 " 

In terms of eigenvalues, we have V(a) = —Nt^cosa so that U(a) = V'(a) = A^Csina. In 
the complex coordinate employed above this corresponds to 

2i \ u 
where u = e ia . 

We search for a single cut solution with A\ = e~ tb and B\ = e tb (See figure 9(a)| ). 



The solution for H(u) is given by ( B.13 ) with the appropriate function h(u). We pause 
to carefully define the function h(u) and enumerate some of its properties. Let u = \u\e ie 
with 6 S (— 7r,7r], then we will define the h{u) as 

h{u) = ({u - e %b ){u - e~ ib )Y . (B.15) 
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The single cut of this function is taken to lie in an arc on the unit circle extending from 
Ai to B\. The function h is discontinuous on this arc with 



ii I o b 9 9 

±2e 2 \ sur sm z - 

2 2 



(B.16) 



where z = e . Away from this cut (i.e. on the gap) on the unit circle 



h{z) = 2ie 2 sin 2 
h(z 



§ - sin 2 § 



2ieT A / s in 2 | - sin 2 § 



>0, 
< 0. 



(B.17) 



(recall that sin 2 | > sin 2 \ everywhere on the gap). On the real axis h(x) is everywhere 
real and is positive for x > 1, but is negative for x < 1. In particular /i(0) = — l. 40 At 
infinity 

h(z) -4 z at |z| — )• oo. (B.19) 
H(u) is now easily obtained from ( B.14 ) as follows. At large z, — — — — 



Nh(z)(z-u) 2iz 



SO 



that the contour integral at infinity in (B.14) evaluates to |. The only singularity in U(z) 
is a pole at zero; using h(0) = — 1 it follows that the residue contribution of this pole to 
(B.14) is S-. Putting these contributions together we conclude that 



H{u) 



c 



in other words, 



$(it) 



1 + 



1 



1 + - I h(u) 
u ' 



u — - 



h(u) 



u 



(B.20) 



(B.21) 



40 In more detail along the real axis 

h(x) = \/(x - cosfc) 2 + sin 2 b > 



x > 1, 



h(x) = -y(a - cosb) 2 +sin 2 b < x < 1. 



(B.18) 
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Using (gjD and flBHj) we conclude that 



ip(a) = — cos — a/ sin sin — 

^ ; vr 2 V 2 2 



(B.22) 



It remains to determine b. At large u, 



h(u) = u — cos b + O 



u 



It follows from ( |B,21| ) that at large 



= £(1 - cosft) + f - J = Csin 2 h - + 0[- 
2 V u I I \u 



Comparing with flB.2|), we conclude 



b A 



sm 



2 C 

In summary the eigenvalue distribution is given by 



(B.23) 



C a A . 2 a 
<Ma) = -cos-W--sm -. 



(B.24) 



We recover the famous GWW [10, 12, 11] result upon setting A = 1 in this result. 



B.2 The capped unitary matrix model 

In this subsection we turn to the study of the capped unitary matrix model. As we have 
explained in the main text, the saddle point equations for this model are given by po + tp 
where po is any continuous function that obeys the boundary conditions and ip obeys the 



integral equation (5.4). ( |5.4|) is a special case of (B.l) with U(a) and A specified in the 
last two of Q5,4| ). Consequently the analysis of the previous subsection applies to this 
case. In particular the eigenvalue density ip is given by flB.4| ) for the analytic function 
= h{u)H(u) where H(u) is given by 



H(u) 



1 



dv- 



U(v) 



(B.25) 



2Wc cuts Nh(v)(v-u) 
(see (lBl3D ) with 

U{z) = V{z) -NV [ % ( W )^, A = 1 - / — P0 (w) (B.26) 
J w z — W J IW 



where both integrals run over the full unit circle (more precisely, the integral can be taken 
to exclude all lower gaps because po vanishes there). 

The contribution to H(u) (|B.25 ) from the second term in U(z) in ( B.26[) (the term 
proportional to V J ^Po(^)f~^) ma y be evaluated as follows. We interchange the order 
of integration in ( B.25| ) to do the v integral first, followed by the w integral. The v integral 
may conveniently be performed by deforming the integration contour to infinity; in that 
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process we pick up the contribution from the poles at v = w. The integral at infinity 
vanishes because the integrand falls of at least as fast as . The contribution of the poles 

at v = w 

J n{w)(w — u) 

where the integral runs over the all cuts and all upper gaps. However it is easily seen that 
the contribution of the cuts to this integral vanishes 41 . In other words the integral in this 
expression receives contributions only from the upper gap, where po(w) = o^r- Putting it 
all together we conclude that 

H(u) = H v (u)+H po (u) 

i f ... n«) 



Hv{u) = ±-J dv 

,J l--- exits 



Nh(v)(v-u) (B.27) 



C it t s 

dw 



HpQ (u) = +i ihtnii ^clc ^ Ju-v L + J_ / 



h(u) ttA J Lugs h(w)(w -u) 

where C cuts is any contour that encloses all n cuts but excludes the point u, L U nit circle runs 
counterclockwise over the unit circle and L ug is the sum over arcs that run counterclockwise 
all upper gap regions on the unit circle. The analytic function $>(z) is given by 

$(n) = + $ P0 (u) 

$ y (V) = h{u)H v (u) = h{u)±- j 



cuts 



Nh{v){v-u) (B.28) 



,, „ , dw , s u + w ih(u) f dw 

<f> po (u) = h{u)H po {u) = iV j — po ( w )— + I 



r , w u — w 7rA It h(w)(w — u) 

^umt circle " ^ugs \ / V J 

The discontinuity in & po is easily evaluated; using ( p.4[ ) we find that across the cuts, 42 

7ih + (7~\ f dm 

*+ (*) - (*) = 4vr^ (z) = —LI / *™ - 4vrp (,). (B.29) 



7T/ 



A JL ugs h(w)(w-z) 



Let us define 

= ^ v (z). (B.30) 

It follows that the eigenvalue density, p(z), of the capped unitary matrix model is given in 
the cuts by 



p(z) = rbyiz) + Vp (^) + Po(z) = Vb v (z) + — -y- j u( ^ f _ _ v (B.31) 



2ih + (z) f dw 

^ X JLugs HW)(W-Z)' 



Of course p(z) vanishes in lower gap regions and equals ^\ m upper gap regions. 



41 The reason for this is that h(w) changes sign across a cut. Consequently the principal value - the 
average of an integral over the upper and lower contour - vanishes over the cut. 

42 Each of the integrals in the third of ( B.27| ) also has a discontinuity across upper gap regions, but these 
discontinuities cancel between the two different integrals in (B.28). As a consequence <& po is analytic across 
all upper gaps. It is also, of course, analytic across all lower gaps. 
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At large u, 



, . ih(u) f dw f dw , , . . 

<^<» = -V / w-w- , - / — PoW + O(-). (B.32) 



7rA h(w)(w — u) J iw u 

It follows that at large u 

^ r(a) = 1 .'mf f +0 (i). (B.33) 

In summary, the saddle point eigenvalue distribution of the capped unitary matrix 
model is given by ( P-31| ) with tpy defined in ( p.30| ) in terms of the function <&{u)y defined 



in (B.28). The contour of integration L ugs runs over all upper gap arcs on the unit circle. 



The locations of these arcs is determined by the normalization condition ( B.33 ). Note that 
the arbitrary function po that we utilized in intermediate steps in our computation has 
disappeared from the final answer as it should. 

Two comments about our solution are in order. First, even though this is not imme- 
diately apparent, it may be verified that the RHS of the eigenvalue distribution function 
( |B.31| ) is real. Second, while the eigenvalue distribution ( B.31 ) is not analytic at branch 
points (boundaries of the cuts) it is continuous. 

2ih + {z) f dw 



nX J L ugs h(w)(w - z) 

vanishes at the boundary of every lower gap (as h(z) vanishes at that boundary). Near the 
boundary of an upper gap, on the other hand, the integral f L h(w)(w-z) diverges, and it 
is possible to demonstrate that 

2ih + (z) f dw 1 



7T/ 



A J L ugs h{w){w — z) 2ir\ 

in this limit. As this is also the value of the eigenvalue density function in the upper gap 
region, it follows that p(z) is everywhere continuous. 



C. Solution of the Capped GWW model 

In this appendix we use the results of Appendix [B] to find an explicit solution to the capped 
GWW matrix model. As reviewed in the introduction to this paper, the capped GWW 
model has four phases. The no gap and lower gap phases of this model are identical to 
that of the original GWW model, reviewed in Appendix ^. The new phases of this model 
are the upper gap phase and the two gap phase. In the rest of this appendix we proceed 



to determine the solution in these phases applying the general formalism of subsection B.2 
above. 
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C.l One upper gap solution of the capped GWW model 

In this subsection we apply the analysis presented above to obtain the solution with a 
single upper gap (but no lower gap) of the capped GWW model, i.e. the special case 



V (u) = — [u - - \ . 

2i \ u ) 



We search for a solution with no lower gap and one upper gap that extends from coun- 
terclockwise from e~ ia to e ia (a is positive). In other words the cut in our problem is the 
complement of the cut in subsubsection |B.1.2| , see Fig. 9(b) . 

The one upper cut solution may be obtained by applying the general method out- 
lined above. The function h(z) = \f (e ia — z)(e~ ia — z) in this section differs from h(z) in 
subsubsection B.1.2] in the location of its branch cut, which is taken to run counterclock- 
wise from e~ M to e ia . We pause to enumerate some of the properties of the function h(z) 
employed in this section. The function we work with obeys 



-hT(e ia ) = h + (e ia ) = + 2ie l 2< sin 



a 



sin 



if 7r > a > a, 



- 2/f ■' -! x / sin 2 — — sin 2 - if 



7T < a < —a, 



hie 1 



-2e l 2 \ I sin 2 - — sin 2 ^ if a > a > —a, 



h(x) = + y (x — cos a) 2 + sin a for x > 1, 



(C.l) 



(x — cos a) 2 + sin 2 a for x < — 1, 
h(z — > oo) — > z, 
h(0) = 1, 

where h~ and h + are the values of function along the branch cut just inside and just outside 
the unit circle. 



Using these properties, the contour integral in the second of ( p. 27 ) is easily evaluated 
and we find 



Hv{u) = Q - 



1 



(1--) 



u 



h{u) 



(the contribution of the contour at infinity is |). It follows that 



h(u)(l--)-(u--) 
u u 



It follows from (jOcb , (|CT|) and dcj) that 



ATTljjy (e 



2 M 



-4C sin — \ I sin" '■ sin* 

s 2V 2 2 



(C.2) 



(C.3) 



(C.4) 



In order to simplify the second term on the RHS of (B.31), simply recall that 

h+(e^) = 2^2, : : - - sin- '- 
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and note that 

It follows from flB.3lQ that 
4vrp(e^) = -4C 



—idw 
h(w) 



dtie** 



2 A /sin 2 § -sin 2 f 



8111 — 

2 



■ 2 M 

sin sm — 

2 2 



2 a\5 2/i+(e^) 



vrA 



d0 cos |(1 -e" 



(sin 2 f - sin 2 f 



sin 2 ^ — sin 2 



-4C 



sm — 
2 



. 2 M - 2 a \ 2 , 2 
Sm 2" Sm 2) + A 



(C.5) 



where the last equality follows from explicitly evaluating the integral. 



We now turn to the evaluation of a. At large u, 3>y(«) = — |(cosa + 1) + 0[M. Now 



( B.33|) at large u gives 
which implies 



C/ \ i 

--(C06O+ 1) = 1 - - 



COS 



2 C 

In summary, the saddle point eigenvalue distribution is given by 



(C.6) 
(C.7) 



p(e*) 
p(e*) 



sm — 
2 



• 2 M ■ 2 a \ 2 , 1 
sm sm — + 

2 2/ 



2vrA 



> a) 



(C.8) 



2vrA 



(M < a) 



where a is defined in ( |C.7| ). 

C.2 Two cut solution of the capped GWW model 

In this section we search for a solution of the capped GWW model 43 that has one upper 
gap centered around a = and one lower gap centered around a = ir. We assume that our 
cuts extend along the unit circle from e m to e lb and from e~ lb to e~ m respectively (a and 
b lie in (0, 7r) with b > a, see also Fig. 9(c) ). In order to obtain the eigenvalue distribution, 
we use the general discussion presented earlier in this subsection, with h(u) as 



h{u) = l(u-e ia )(u 



-ib 



(C.9) 



h(u) is take to have cuts along the arcs enumerated above. We pause to summarize some 
of the properties of the analytic function h(u). 

Along the unit circle outside the cuts, the h is analytic as 



h(( ! ") Ir'" v / (sin 2 | - sin 2 ^)(sin 2 h - - sin 2 



-a < < a, 



hie 1 



Ae l \l (sin 2 - - sin 2 -)(sin 2 - - sin 2 - ), 



(CIO) 

-7T < 9 < ~b, b < 9 < 7T. (C.ll) 



* As in the previous subsubsection V'(a) = (^Nsina. 
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Along the cuts, h becomes discontinuous with 



± Me %t \\ (sin 2 ^ - sin 2 ^)(sin 2 \ - sin 2 a < 9 < b, 



T Ue ie \l (sin 2 °- - sin 2 |)(sin 2 ^ - sin 2 ^ 



-b < 9 < -a, 



(C.12) 



(C.13) 



where h = h + on the outside of the circle and h = h on the inside. Moreover h is real 
and positive along the real axis; explicitly for real x 

h{x) =\J (1 — 2x cos a + x 2 )(l — 2x cos b + x 2 ), — oo < x < oo. (C14) 



Note in particular that h(0) = 1. At large u we have 

h(u) ~ +u 2 — u(cos a + cos 6), u ~ oo. 
Proceeding as in the previous subsubsection we find 



(C.15) 



H v (u) 



1 



u 



and 



which follows that 



h(u)C 



1 



u — fr 



/i(n) 



4mpv (u) 



-h + {u)- 
u 



The full eigenvalue density can now be computed using (B.31). The second term on the 
RHS of ( B.31| ) can be simplified by moving back to angle variables. Substituting w = 
e td , u = e ia , converting the integral over w to integral over 9 and using ( p,17| ) below we 
find the following explicitly real expression for the density distribution 



47rp(a) 
h 



l sin a l 

7rA 

a 



sin sin — 

2 2 



• 2 ■ 2 a \ i 

sin sin — ii 

2 2 ' 

d0 



where 



(C.16) 



sin" I - sin 2 |) (sin 2 | - sin 2 £ 



'-" (cosW - cosn ) J jsin 2 £ - sin 2 f) (s 
It remains to evaluate a and b. Setting the coefficients of terms of 0{u) and 0(1) 



in ( B.33 ) to zero (and substituting w 
equations 



1 



4ttA 
1 

4^A 



da 



1 



,2 a 



da 



cos a 



e ia to simplify the integrals) 44 we obtain two 

c 



,2 a 



sin 2 | - sin 2 § w sin 2 | 



sin 2 | 



1+2 ( cos a + cos ^) 



(C.17) 



which determine a and b. ( |C.16 ) and (C.17) provide a complete solution of the two cut 
model. 



Under this substitution h(w) — 4e la wsin 



- 2 | - sin 2 f a/sij)-' 



sin 2 # where the square root is positive. 
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C.3 Special Limits of the capped GWW model 
C.3.1 p{a) as b — > ir 

Let us first consider the eigenvalue density in b — > it limit. Setting b = it in ( p,16| ) and we 
obtain 

Ai:p(a) 



sin a cos ■ 



2 a o a 

sin sin — 

2 2 



d.6 



a (cos I 



cos a ) cos [sin 2 | — sin 2 |^ 



(C.18) 



Using 



f 



dO 



7T 



(cos 9 — cos a) cos f \ (sin 



,2 a 



,2 



9 (X 
COS^ ^ 



sin ^ w sm 2 2 — sin" | 



2 a COS 2 



together with (7.7) we find 

p(a) 



1 



c 



2vrA vr 



a 



sm ■ 



2 a 2 a 

sin sm — 

2 2 



(C.19) 



(where C and a are given in (fr77|)) in perfect agreement with Q7.3|) . 
C.3. 2 p(a) asa->0 

Now let us examine the a — > limit of ( |7.6|) . While the limits of integration coincide in 
this limit, and the integrand diverges in a precisely compensating manner, leading to a 
finite result. The integral is easily evaluated in this limit; combining the result with the 
first equation in (|7.5| ) we find 



C _„ a . L 2 b ■ 2 a 



p(oi) = — cos — \j sin sin — 

rv ; vr 2 V 2 2 



(C.20) 



(where £ and A satisfy (7.9)), in perfect agreement with (\l.2\). 



C.3. 3 The large C limit 

In this section we elaborate on the evaluation of the integrals appearing in Eq.( C.17| ),( |C.3C| ). 



To start with, let us consider 



1 = 2 1 d9- 

'° */(sin 2 (§) 



sm 



2 (D) 



,2/6n 



sm 



2/6* 



D) 



(C.21) 



Since b is very near a in the large £ limit, the integral receives most of the contribution 
from the integration region near a, and it diverges. To separate out the divergent piece, 
let us divide the integral into two parts I = I\ + 12 , 

1 



h=2 I dd 

'"-7 ^(sin 2 (f ) - sin 2 (f )) (sin 2 (|) - sin 2 (f )) 

ra-7 1 



h 



(C.22) 



(10 







sin 2 (f ) - sin 2 (f )) (sin 2 (|) - sin 2 (f )) 



- 69 - 



where I2 is finite whereas 1\ contains a divergent piece. Our aim is to compute a, b in terms 
of C- Now we take b = a + -X, where M is large. We make a change of variable 9 = a — %ry, 
and expanding Ji and I2 in O (6> 2 ) and O (57) 2 respectively we obtain 

h = 8csc(a)sinh _1 (y/ltj + j- (2 y/M (M + 1) cot a esc a - 4 cot a esc a sinh" 1 \/m) 
+ O ( 7 2 ) 

^2 = 4 esc a ^log sin ( a — tt) — 1°S s i n (1^ 

7 / sin(a — 7) / /7\ / 7 

+ — esc a 2 7 r + 2 cot a log sin — — log sin a — — 

M \ cos a — cos(a — 7) V V2/ V 2 

+ °(f) 2 - < C23 > 

When re-expressed as a function of e(= -X) and M, /1 + I2 must (of course) be independent 
of M. Within the perturbative expansion presented above this works as follows. Terms of 
order 0{^j) in the expansion of /21 when re-expressed as a function of M and e, yields 
expressions of order -i, and e 2 (along with logarithmic corrections). As we have not 
computed these second order corrections, it follows that the leading M dependence of the 
O(l) part of our answer must be at order j^, and that the leading part M dependence of 
the O(e) part of this integral occurs at O(jj); we have explicitly checked that this is the 
case. These sub leading M dependences in the answer are an artifact of our perturbative 
expansion, and would vanish if we carried out our expansion to higher orders. The correct 
answer for the integral upto 0(e) is simply given by taking M to infinity in the coefficients 
of the O(l) and 0(e) parts of our integral. 45 We finally arrive at 



I = h + l 



2 



sin 2 (a) V ° V 6 J J sin(a) 

In a similar way one perform second integral appearing in Eq.( |C.17[ ) we obtain 

cos(6>) 



ra 

2 / dO 



„ / ( a - m 2 (%) - sin 2 ^^ fsin 2 ^ _ a ^2, 



(§)) (sin 2 (|)-sin 2 (f)) 

/ , /8sin(a)\ 1 / , /8sin(a)\\\ ^ „. 

= 2 (2 a + 2 cot(a)log {-11 j - (-1 + log {-11 j j j + 0[«*|. 

(C.25) 

Now using J,/3 in ( |C . 1 7| ) we get two equations for a and e. Solving these equations we 



45 Note that taking M to infinity in the coefficients of higher order terms - e.g. the coefficient of e 4 - 
would yield a divergent result at this order in perturbation theory. This is because we have not yet included 
the third order terms in I2 that would cancel this divergence. As a consequence, the correct procedure to 
estimate the integral expanded to first order is to truncate I\ + 12 to first order in e and then to take M to 
00 in the final answer. 
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obtain 

e = 8sin(vrA)e- m:Asin( " A) + 0[< e - 2 < Asin (- A )] 

a = 7rA - 4sin(^A)e^ Asin ^ A ) + 0[( e~ 2 < Asin ^ A )] / Q 2g ) 

= *A-|e + 0[C e 2 }. 

Now in a similar way, one can compute the leading term in the eigenvalue density using 
Eq. flC.ldf ), In order to proceed we make following change of variables, 

9 = a — ey, b = a + e, a = a + ea\. (C.27) 

Note that a\ takes value between to 1. The value of p takes the form 

4vrp = — - cos -1 Ja{. (C.28) 

As a check of our result, we show J pda = 1. The integral of p can be written as 

/7r rb r-a 

pda = 2 / pda + / c?apo 
7r J a J —a 

= 2 / epdai + / ( J,/ a ^A (C29) 
+ (1 



2?rA v 2vrA y 
= 1. 

C.3.4 Eigenvalue density near the end points of cuts 

In this subsubsection we study the behaviour of the two cut eigenvalue density function 
p(a) given in ( p,16| ) as a tends to the edges a and b of the distribution function. We work 
at fixed C and A (hence fixed a and fixed b). 

Let us first consider the limit a — > b. Clearly, to leading order 



4vrp(a) = B^Jb-a + 0((b - a 



I Sm |l // ■ 2 b . o a \ ■ h b 

B = -^ry {sm 2" Sm 2 )sm 2 C ° S 2 (C.30) 



dd 



a Z{sm 5 - sin 2Ja/v S 2 _ 2-'V sm 2 _ 2> 



2 uxx± 2 / \ / xx± 2 o±±j. 2 / 2 2 / 

It is not difficult to check that the coefficient of \/b — a agrees with the coefficient of the 
same term in the Gross- Witten-Wadia solution when a is taken to zero, and that it vanishes 
when b — > ir. 

We now study the limit a — > a. The estimation of the eigenvalue density function in 
this limit is more subtle; we follow the method described in subsection C.3.3| . Let a — a = e. 
The expression for p(a) in ( |C ■ 16| ) may be recast as 

p(a) = h + h 
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where 



sin(a)^/ (cos(a) — cos(b)) sin(a) 
I\ = : 



7TA 



d6- 



a-7 (cos(0) - cos(a)) A /(sin 2 (f ) - sin 2 (f )) (sin 2 (f ) - sin 2 (|)) 

' " ' ' ' (C.31) 



sin(a)-v/ (cos(a) — cos(b)) sin(a) 1 
h = t e 2 

a— 7 

d0 



(cos(x) - cos (a) )^ (sin 2 (f ) - sin 2 (|)) (sin 2 (|) - sin 2 (f )) 



As in the previous section we have divided the integral in ( C.16| ) (which can be taken to 



run over the range (0, a) into an integral from (0, a — 7) and an integral from (a — 7, a)). 
We take 7 = Me where M is a large but number that is held fixed in the limit a — > a i.e. 
e -»• 0. 

Inside the integral in I\ , it is legitimate to Taylor expand a and the integration variable 
9 about a. The contributions of successive terms in this Taylor expansion are suppressed 
compared to leading terms by increasing powers of e. It is not difficult to verify that 

h = (~ + hi(M) \ + e h 2 (M) + 0(e§) 
fci(M) = -^= + 0(-^-) (C.32) 



fi2(M) = — — rr (1 + 5cos(2a) — 6 cos a cos 6) . 

27rA(cos a — cos b) sin a 

The expression for I2 may be processed as follows. Of course I2 = (h — h) + h, where we 
choose I3 as 



sin(a)-v/(cos(a) — cosT&Y) sin(a) 1 n a ^ 4 

I 3 = v — —ei I dx. 

ttA Jo y'cos(a) — cos(fo)-y/(cos(a) — cos(x)) 3 

(C.33) 

The point of this manipulation is as follows. In the expression for I2 — I3 we first subtract 
the integrands before performing the integral. In order to proceed further we proceed as 
follows. 

h-h 

^sin(a) 1/ (cos(a) — cos (6)) sin(a) 



'e- 



7rA 



a-7 f 4 4 \ ! (C.34) 
H dx 

y ycos(a) — cos(6) y'cos(x) — cos(6) y y / (cos(a) — cos(x)) 3 

1- . rrr sin a _ 3, 

7rA(cos a — cos 0) 
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where 

1 ( a ) \J (cos(a) — cos(fo)) sin(a) 



A Sml 

A = 



VTA 

4 4 ^ 1 

H eto. 



(C.35) 



V y cos(a) — cos(6) y / cos(x) — cos(6) / (cos(a) — cos(x)) 3 
The I3 integral can be performed easily and is given by 

/ 3 = -/ii(M) + h 3 (M) - — \' l2;; + 0[e5] 



vrA^cot (f ) (c>36) 



, /rr3cota ^ r 1 
/i 3 (M) = -eVM — + Q[ 



\ ' - L 3 J' 

where symbol E refers to elliptic function of second kind. Adding all up, we finally obtain 
eigenvalue distribution as a — > a to be 

1 1 3 

p(a) = — - H D\/a — a + 0(a — a) 2 

27TA 4ir 

with D _ A _«(«i^rt» < c - 37 > 



We notice once again that near a, the eigenvalue approaches -^j as ea. As a check of our 
result, we take b = ir. In this case result should reduce to that obtained from one cut 
solution. In this limit the A in ( C.35| ) can be easily computed and is given by 



4£(a|csc 2 (§)) (sin a) 

J\ — 



vrAJcot (f ) v^Acos 2 ! 



AE (a I esc 2 (% 1 j </ . a a 



-AC sin — \ I sin — cos — , 
s 2 V 2 2' 



(C.38) 



cot (§) 

where in the last line we have used ( |7.8| ). So finally we obtain as 



a 



1 C a a a 3 
P(a) = 2^A " ^ Sm 2 V Sm 2 C ° S 2 + ° {a ~ a) ~ 2 

which is in perfect agreement with ( C.19| ) as we take a — > a. 

C.4 Level rank duality of the solution to the capped GWW model 

In this subsection we will directly verify that the explicit solution of the capped GWW 
model obeys ( |7.13| ). We will find the following notation useful. We denote the 'no gap' 
saddle point eigenvalue distribution of the capped GWW model by /3 n5 (A, £, a). In a similar 
manner we denote the one lower gap, one upper gap and two gap eigenvalue distributions 
by pi g (X,C,a), p„ 9 (A,C,a) and p tg (X,C,a) respectively. 
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No Gap Solutions: In the no gap phase 



Pngi^C, 



n ) 



1 + C cos a 
2n ' 



Using flagp it is easy to directly verify that 



Pn 9 ((1 - A), a) = ^A- (J- _ pn9 (A, C, a + tt)) 



1 - A V2ttA 

so that level rank duality maps the no gap phase to itself. 
One gap solutions: In the one lower gap phase (see ( |7.2| 



(C.39) 



(C.40) 



Plg(\C,a) 



Ccos(f) fl . 2 a . 2 a 1 

- ' - — sur — for sin — < - 



TT 



c 



c ■ 2 a 1 

for sin — > 



c 

I 

C' 



(C.41) 



On the other hand in the upper gap phase (see ( |7,3|) ) 

m.rr- 



1 

2^A 
1 



sin 



7T 



c 



COS 



9 a 2 a 1 — A 
z — tor cos — < 



for cos 



> 



AC 
1 - A 



(C.42) 



2ttA 2 AC 

Using ( C.42| ) and ( p.41| ) it is straightforward to directly verify that 



Pig ((1 - A), y-^- 



CA 



Pug [{I ~ A), j— j 



a 



a 



A 



1 - A V 2vrA 



1 



A 



1 



Pug (A, C, a + 7T) 
ft 9 (A,C,a + < 



(C.43) 



1 - A V2?rA 

It follows that level rank duality exchanges the one lower gap and one upper gap solutions. 



Two gap solution: The two gap solution is listed in (7.6) with the end points a and b 
given in ([7.5D . We believe it is true that 

p tg ((1 - A), a) = A_ f _L _ ptg (A, C, a + tt)) , (C.44) 

so that level rank duality maps the two gap solution to itself. 

We have not succeeded in directly verifying ( |C.44j ) starting with ( |7.6j ). However we 
have generated some numerical evidence that ( |7.6| ) obeys ( |C.44| ), as we now proceed to 
describe. 

a and b, the end points of the gaps in the two cut solution, are both functions of A and 
C- As we have explained in Appendix C.5 it is formally useful to invert this dependence 
and regard A and C as a function of a and b. In particular the dependence of A(a, b) on a 
and b is listed in the first of ( C.5C ). Level rank duality interchanges upper and lower gaps 
and so requires that 

A(vr-6,vr-a) = 1- A(a,6). (C.45) 
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Using the explicit integral representation for A(o, b), ( |C.50 ) we have numerically verified 
( C.45| ) to high precision for several randomly selected values of a and b. 

We have also generated numerical evidence for the self duality of the two gap eigen- 
value distributions themselves in the following manner. Our solution for the eigenvalue 
distribution (7.6) takes the following structural form 

2TT\p tg (\,(,a) = f d0v(a,b,a,9), (C.46) 

J —a 

where the explicit form of the function v may be read off from (|7,6j). On the RHS of 
( C.46j ) a and b are the functions of A and ( determined by ( |7.5|) . As we have explained 
above, however, for the purposes of the current section it is more convenient to invert this 
dependence and regard A and £ as functions of a and b. Let us define 

H(a, b, a) = 2ir\(a, b)pt g (X(a, b), ((a, b),a). 

( C,46j ) may be rewritten as 

H(a,b,a)= I d0v(a,b,a,9). (C.47) 

J —a 

Now the prediction ( |C,44j ) may be rewritten as 

H(a, b, a) + H(tt - b, n - a, a - ir) = 1. (C.48) 

Using the explicit integral representation for ( p.47 ) we have numerically verified ( |C.48D to 
high precision for several randomly selected values of a, b and a. 

While the numerical evidence for ( |C.45 ) and (C.48) is impressive, it would certainly 
be useful to find an analytic proof of these equations. We leave this for future work. 

Level Rank duality of the phase transition points: As we have explained above, 
for A < 2> capped GWW model undergoes two phase transitions at 

Cf (A) = 1, 

for no gap to lower gap and 

for lower gap to two gap respectively. On the other hand, for A > |, the two phase 
transitions occur at 

tfw - V- 

for no gap to upper gap and 

c£(A) 1 



4A(1 - A)' 

for upper gap to two gap respectively. Under level rank duality a value of A < ^ maps to 
a value of A greater than i. Level rank duality requires that 

A C f(A) = (l-A)Cf(l-A), 
A<|(A) = (1 - AKtU - A). 

( C.49| ) is very easily verified using the explicit expression for Q, tf, and presented 
above. 
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C.5 Behavior of (A, £) with respect to (a, b) in the one lower gap and one upper 
gap solution 

The set of equations ( C.17| ) is equivalent to the following set 

A( a ,6) = — / da (G(a,b,a) + 1 — - ) , 
4vry_ a V G{a,b,a)J 

1 r 1 (C50) 

A(a,6)C(a,6) = — / cfa- 



where 



'" • " /-n 2 f - sin 2 f Jsin 2 § - sin 2 § 



sin 2 | — sin 2 ^ 



G(a, b, a) = v =. (C.51) 



sm 2 ^ — sin 2 ^ 



One can confirm it by subtracting ^ (cos a + cos 6) times first line from the second line in 
( C.17| ). First let us observe the behavior of A with respect to the change of b with fixed 
a = a Q . In the case ir > b > b' > a Q > 0, since 1 > sin 2 | > sin 2 ^ > sin 2 ^ > 0, it becomes 

G(a Q , b, a) > G(a Q , b', a) > 1. (C.52) 

Therefore 

G(a Q , b, a) + — \ ) > (G(a , b' , a) + — 1 ) > 2. (C.53) 
&(a ,b,a) J \ G(a ,b',a)J 

Hence in the fixed a , A(a , b) is a increasing function with respect to b, i.e, 

A(a G , b) > A(a D , b') at vr > b > b' > a Q > 0. (C.54) 
The range of A in the fixed a = a Q is 

— <A(a ,k)<l--cos— (C.55) 

7T 2 2 

where 

A(a ,a ) = — , A(a ,7r) = 1 - - cos -A (C.56) 

7T 11 



From the factor 1/ysin 2 | — sin 2 ^ in the second equation of ( C,50| ), we can immediately 
see that AC(a , b) is the decreasing function with respect to b. Moreover, since the A(a , b) 
is increasing function, ((a ,b) becomes decreasing function. The range of ((a ,b) is 

w - ¥ L__< aao , b) < co , (C.57) 

where ^ 

CK, vr) = —- j — — , C(a Q , a ) = oo. (C.58) 

2 (1 - i cos -f) cos -f 

We can also see that 

««., > ««., «) > min ( 2(l _^ f)cos J = I- (C59) 
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Figure 10: 



Plotting 



(a) 

(AM), CM)) 



with the values of (a, b) 



(f'|)'(f'f")'(f'f")'(f'¥) and (f'¥)- From this § ra P h ' we can S uess that are in " 

creasing functions of a, and we can also see that A is surely a increasing function of b and £ is a 

decreasing function with respect to b. The region £ > 4J2, A < ^ and £ > 4A( - 1 1 _ A - ) , A > \ is shaded. 

All points are located on the inside of the shaded region. 



Hence £ must be larger than 1, and we have exactly shown that £ has a unique minimum at 
a = 0, b = ir. As a summary, A is a increasing function of b while £ is a decreasing function 
of b. We can also guess the behavior of A and £ with respect to a, with fixed b = b Q by 
a numerical calculation. (See Fig. 10(a) ). From this, we can guess that A is increasing 



function of a and the also £ is increasing function of a. With fixed a, let us consider the 
one parameter line (A(a D , b), £(a Q , b)) with parameter b in the A — £ plane. (See Fig. 11 (a)] ) . 
Since the X(a ,b) is increasing function with respect to b while £(a G ,6) is the decreasing 
function, £ = £(A) behaves as the decreasing function with respect to A. The £ is divergent 
at 

£(a ,6)-»oo at A(a G ,6)->— (6— > a Q ). (C.60) 

7T 

In the fixed b = b Q case, we can also consider the one parameter line (A(a, b ), £(a, b a )) with 
parameter a in the A — £ plane. In the case of fixed b = b = it, the one parameter line 
(A(a, 7r), £(a, tt)) is represented as 



c 



1 



4A(1-A) 

On the other hand, in the case of fixed a = a 
is represented as the 

^ 4A 2 



I<A<1. 



(C.61) 

0, the one parameter line (A(0, 6), £(0, b)) 



0<A<I. 



(C.62) 



Next let us consider the line (A(a G , b), £(a Q , b)) with fixed a Q / 0. In this case, due to 



C(a ,7r) 



1 



4A(a ,7r) (1 - X(a ,7r)) 
((a ,b) > £(a ,vr) = 



and h < X(a ,b) < X(a ,ir), we can see that 



> 



4A(a G , tt)(1 - A(a , tt)) 4A(a D , 6)(1 - X(a , b)) 



> 1. 



(C.63) 
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<T=1/(4A(1-A)) 

(A(a, 7i), ((a,n)) 



(a) 

Figure 11: The schematically drawn two lines (A(a Q , b), ((a , b)) and (X(a' , b), ((a' , b)) with fixed 
a a and a' a with a < a' . They are drawn by dashed line. For each line, £ becomes a decreasing 
function of A. 



for b < tt. Then the lines (A(a , b), C(a , b)) are located on the inside of the region £ > 
4A(i-A) wnen -^( a o,b) > 1/2. Also for the case \(a ,b) < 1/2, we can show that the 
lines (A(a Q , 6), £(a , 6)) are located on the inside of the region ( > 4j. To show it, we 
should note that the value of (a, b) is uniquely determined if we specify the 2 parameter 
(A, £). This means that there are not any intersection between two lines, (A(a D , 6), £(a , b)) 
and (A(a' G , 6), C( a ' i b)) with fixed a Q ^ a' Q . In particular, there is no intersection of line 
(\(a ,b),C(a ,b)) and the (A(0, 6), C(0, b)), which is written as ( = For the lines 

(\(a ,b),C(a ,b)) with fixed < a < |, it becomes 

C(a ,b) > 1= 1 at A(a ,6) = i. (C.64) 

4A(a o ,0) / 2 

So for A(a G , 6) < |, the line (A(a D , 6), £(a , &)) must be located inside the £ > jj^ since the 
line can not have the intersection with the line (A(0, b), £(0, &)). 46 This shows that two 
gap solution exists if and only if £ > and £ > • ^* ur numerical result also shows 



it. (See Fig. p)(al) 



D. Level-rank duality of the saddle point equation in the multi trace po- 
tential 

We consider the multi-trace potential, 

/ N \ m 

y = ^^4,n(W n ) m + c.c = ^^i m ,„ +c.c. (D.l) 

m n m n \i=l / 



40 



For the line (A(o , b), C(flo, &)) with a > -| we do not have to worry about the region with A < ~. 
Because the range of the A of the line is 1/2 < a /n < A < 1. 
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The corresponding potential in the dual theory is obtained by tiU n o (— 1)™ +I tri7 n , 

/k-N \ m 

y = ^^A m>n (tr(-ir +1 C/"r + C . C = ^^^ m , n [£- e ***+'0j +cc (D.2) 

m n m n \ i=l / 

Derivative of the potentials with respect to an eigenvalue a/ are 

/ N \ m ~ l 

V ( = da~ V = E E E ^ {imne ma >) + c.c, (D.3) 

1 m n \i=l / 

rl /k-N \ 771-1 

V{ = j-V = E Am > n ( E -e in(Ql+ " } ) (-imne m(ai+l) ) + c.c. (D.4) 



i=l 



In ( p. 3D , please note that only the factor (imne mai ^ carries the label I dependence of the 

/ . \m— 1 

V l (there is no longer label I dependence in the I YLi=i e mai I because it is the sum over 



the label). Similarly, in (D.4), the label I dependence in V{ is carried only by the factor 
(-jmne m ' ai+?r '). Saddle point equations are obtained for each eigenvalue ai by taking the 
derivative as, 

° - £ ("-£- h^)) - * - 5> (^) • < M > 



da>i 



V - J> (2*. 2^)) = «? - E «* (^r 5 ) • < D6 > 



Now let us consider the large N limit in the original theory. When we take the large N 
limit, we first fix the residual Weyl permutation, exchanging the order of the eigenvalues. 
After fixing the Weyl permutation, the order of the eigenvalues obey following 

< a\ < a 2 < . . . < o/v < 2ir. (D.7) 

Then there is a one to one relationship between the value of the eigenvalue a\ and the label 
I such that 

i < I O at < ai. (D.8) 

In the large N limit, these labels can be regarded as the continuum variable x = -L, (I = 
I,... N). Moreover due to the flEXf) , /3(x) = a u (x = jj) becomes the increasing continuum 
function with respect to x, so 

I i 

I < i45 — = x < x' = — <^ ai= P{x) < f}(x') = ai. (D.9) 
And from these facts, the trace is written by the integration over the (3(x) as, 

V N / dx = N / — d/3 = N / dPptf). (D.10) 
i=1 Jo Jo a P Jo 
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From (|D.9| ), we would like to emphasize that value of eigenvalue variable f3 represents the 
label I of the eigenvalue a\ = j3 in the discretized representation. According to these, V'((3o) 
is just a large N representation of the V( with /3q = P(jt) = &i, 



V ( = E E A ^n E ^ {imne 1 ™*) + c.c 

m n \i=l J 

-»• *"(A>) = E E / dPp{f}y^ J (imne^J + c.c. (D.ll) 

So, with keeping in mind that the variable j3 represents the label of eigenvalue, V'(/3o + ir) 
becomes the Vl with the label V such that ay = a\ + tt = /3q + 7r, 

/AT X™- 1 

V'(A> + tt) = V{, = E A ^ E ema " (imne iB «" ) + c.c 

m n \fc=l / 

/ f-2-K \ m-1 

= £ j4m '« ( N J d(3p((3)e m n (imne in V°+^) + c.c. 



m n 



(D.12) 

This is also the case in the dual theory. At the large k — N limit, V'(/3o) is 



V'(po) = ^5> m ,„ E {-imne ma ')+c.c 

m n \ i=l / 

( r 2iT \ m ~ 1 / \ 

= E Am > n [-( k - N ) / d(3p((3)e in ^ +n ^ J (-imne in ^ 0+7r n + c.c, 

(D.13) 

with /3o = ct[. 

Now let us return to the large N limit of the saddle point equation and the discussion 
of the level-rank duality. The large N (or large k — N) representation of the saddle point 
equations are 

y l = E cot (^ir 1 ) ^ v ' m = NV S dpm cot i^Y 1 ) ' (D - 14) 

V> = £ cot (2^) -+ V'(Po) = (k - N)V J dPpifi) {^Y 1 ) ■ (D-15) 



As we did from ( |5.1C| ) to ( |5. 13 ) in the subsection ^2, assuming that p(/3) is the solution of 
the fling ), following j5(0) 



k 2ttA 

is a solution of the (D.15) for the dual theory if 



A 1 ' p(J3 + *)) (D.16) 



V'{p ) = -V'(Po + ir) (D.17) 
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is satisfied with ( p,16|) . Then we only have to show that (|D,17| ) is satisfied with (|D,16| ). 
We can check the statement ( p.17 ) directly. From ( D.13| ), by substituting ( D.16| ), we can 



check 

(■2-k \ m ~! 



V'{J3) = Y<Y1 Am > n ~ N ) df3p(f3)e in(l3+n A (-imne*^ 4 ^) + c.c. 

m n ^ J® ' 

( f 27T \ m-1 

= An,™ f AT / d/3p(f3 + 7r)e in ^ +7r ) J + c.c. 

= - 22 ^m,n C AT ^ d/3p(/3)e in/3N ) (imne ! "(' 3o+ ^) + c.c. (D.18) 



Comparing ( p. 18 ) with the ( D.12| ), we can see (D.17) is satisfied with ( D.16| ). 



Then this concludes the proof of the level-rank duality of the saddle point equation in 
the multi-trace potential terms. 

E. High temperature limit of the partition function of a gas of non renor- 
malized multitrace operators 

In this appendix we compute the partition function of a gas of non renormalized multi- 
trace operators, in the high temperature limit, for the theory of regular bosons, regular 
fermions and the super symmetric theory with a single chiral multiplet in the fundamental 
representation. 

The single letter partition function for bosonic theory is given by 

I 1 + x 
(1 - xf 



(E.2) 



where x = e r . In the large temperature limit we obtain 

x = l-^ + oQ^j, f b (x) = 2T 2 + O(T ). 
The multi-trace partition function is given by 

oo 1 

In Z = Y j -f b (x n ) 2 (E.3) 

71=1 

which in the large temperature limit is given by 

oo 1 

In Z = 4 T 4 2 — + 0{T 2 ) = 4 T 4 C(5) + C(T 2 ). (E.4) 



71=1 



InZ in (|E.4| ) gives the value of partition function for the gas of multi-trace operators of 
free bosons on the sphere of volume V2 = 47r. 

For the fermionic theory, the single letter partition function is given by 
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where x = e t . In the large temperature limit we obtain 

f f (x) = 2T 2 + O{T ). 
The multi-trace partition function is given by 

oo 1 



n=l 



(E.6) 



(E.7) 



which in the large temperature limit is given by 

OO j 

lnZ = 4T 4 ^ — + 0(T 2 ) = 4 T 4 C(5) + C(T 2 ), 



(E. 



71=1 



which is same as that of boson obtained in (EL4). InZ in ( E.8 ) gives the value of partition 
function for the gas of multitrace operators of free fermions on the sphere of volume V2 = 47r. 

For the M = 2 supersymmetric theory with a chiral multiplet, the multitrace partition 
function is given by 

00 

lnZ = E-(/^ n ) + (- 1 ) n+1 //( 2;ri )) 2 



71=1 



which in the large temperature limit is given by 

00 00 1 01 

lnZ = 8T4 E-T + 8T4 E(" 1 ) n+1 ^5 = - T 'C(5) + 0(T* 



71=1 



71=1 



71" 



(E.9) 



(E.10) 



InZ in (E.10) gives the value of partition function for free M = 2 supersymmetric theory 
with a chiral multiplet on the sphere of volume V2 = 4-7T. This exactly matches with first 
term in (6.21). 
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